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ABSTRACT 


Bistatic  radar  has  some  properties  that  are  distinetly  different  from  monostatie 
radar.  Reeently  bistatie  radar  has  reeeived  attention  for  its  potential  to  deteet  stealth 
targets  due  to  enhanced  target  forward  seatter.  Furthermore,  the  feasibility  of  hitchhiker 
radar  has  been  demonstrated,  whieh  allows  passive  radar  reeeivers  to  deteet  and  track 
targets.  This  thesis  developed  a  software  simulation  paekage  in  MATLAB  that  provides 
a  eonvenient  tool  to  examine  the  bistatie  radar  design  parameters  and  prediet  system 
performanee.  The  software  model  is  suitable  for  instruetional  purposes  due  to  its  user- 
friendly  graphieal  user  interfaee.  Several  bistatie  radar  applieations  were  used  to 
illustrate  the  software  features,  and  their  results  were  analyzed  and  diseussed. 
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EXECUTIVE  SUMMARY 


Bistatic  radar  uses  different  antennas  separated  at  a  substantial  distanee,  for 
transmission  and  reeeption.  This  fundamental  differenee  in  the  bistatie  radar  has 
generated  several  properties  that  are  distinetly  different  from  the  monostatie  radar.  These 
unique  properties  have  offered  bistatic  radar  certain  advantages  in  several  eommereial 
and  military  applieations. 

Reeently  bistatie  radar  has  reeeived  attention  for  its  potential  to  deteet  stealth 
targets  due  to  enhaneed  target  forward  seatter.  Furthermore,  the  feasibility  of  hitehhiker 
radar  has  been  demonstrated,  whieh  allows  passive  radar  reeeivers  to  deteet  and  traek 
targets,  using  transmitter  of  opportunity.  Thus,  it  is  likely  that  bistatie  radar  will  play  a 
much  bigger  role  in  near- future  military  warfare. 

Design,  development,  and  testing  of  radar  systems  have  traditionally  been 
expensive  and  time-eonsuming  tasks.  Radar  designers  are  often  faeed  with  many  design 
iterations  to  optimize  the  system  performance.  Thus,  using  eomputer  simulation  and 
modeling  will  provide  the  designer  the  essential  tools  to  examine  design  parameters  and 
to  evaluate  the  system  performance  before  the  aetual  implementation  in  hardware. 

There  are  many  software  applieations  available  for  monostatie  radar  simulation 
but  none  exits  in  the  publie  domain  for  bistatie  radar.  The  objective  of  this  thesis  is  to 
develop  a  general  software  model  to  analyze  the  bistatie  radar  system  performance  as  a 
funetion  of  the  system  design  parameters. 

The  thesis  foeused  on  the  design  and  development  of  a  Graphieal  User  Interfaee 
(GUIs)  in  MATLAB  to  integrate  the  bistatie  radar  parameters  and  to  examine  several 
eritieal  performanee  tradeoffs.  The  end  produet  is  a  program  named  Bidarsa,  whieh 
provides  a  eonvenient  instruetional  tool  for  engineering  eommunities  to  understand, 
analyze,  and  predict  bistatie  radar  performance. 

The  thesis  presents  an  overview  of  bistatie  radar  definitions  and  the  parameters 
involved  in  the  range  equation.  The  radar  deteetion  eontour,  the  ealeulation  to  determine 
target  loeation,  eonsiderations  for  the  bistatie  footprint,  the  doppler  relationship  for  a 


XV 


moving  target,  and  the  eleetronic  countermeasures  (ECM)  issues  are  discussed  with  their 
supporting  mathematical  equations.  These  equations  were  used  to  formulate  the  software 
model. 

The  Bidarsa  program  consists  of  eight  major  modules  developed  and  integrated  to 
carry  out  the  required  performance  prediction  and  system  analyses.  A  set  of  radar 
parameters  was  compiled  to  demonstrate  the  software  functions.  The  thesis  further 
illustrates  the  software  functions  with  two  bistatic  radar  applications.  The  first 
application  used  bistatic  radar  to  monitor  shipping  channel,  and  the  second  application 
investigated  the  feasibility  of  using  direct  television  broadcast  satellites  and  ground-based 
bistatic  receivers  to  detect  targets. 


I.  INTRODUCTION 


Bistatic  radar  has  come  a  long  way  from  its  initial  days  well  before  the  outbreak 
of  World  War  II.  Many  seientists  were  trying  new  radar  configurations  where  the 
reeeiver  was  located  separately  from  the  transmitter.  The  main  motivation  at  that  time 
was  to  gain  sufficient  radio  frequency  isolation  between  the  two  electronics  sets, 
partieularly  with  the  numerous  ongoing  eontinuous  wave  (CW)  experiments  [1].  This 
fundamental  differenee  in  the  bistatie  radar  over  monostatic  radar  has  offered  many 
opportunities  in  both  commercial  and  military  applications.  Over  the  years  and  through 
three  resurgences,  many  bistatic  radar  systems  have  been  proposed,  developed,  and  tested 
in  many  eountries  including  the  United  States,  United  Kingdom,  France,  Soviet  Union, 
Germany  and  Japan. 

This  unique  difference  has  offered  bistatic  radar  certain  advantages  for  particular 
tasks  [2]; 

•  The  receiver  is  completely  passive,  and  hence  undeteetable,  by  the  electronic  support 
measures,  and  it  is  safe  from  attaek  by  anti-radiation  missiles  or  deliberate  directional 
interference  and  jamming. 

•  The  reeeiver  ean  be  located  in  a  favorable  area  where  transmitters  are  not  allowed, 
sueh  as  flammable-liquids  stores,  gas  terminals,  etc. 

•  No  transmit-receive  switch  or  duplexer  is  required.  These  devices  are  lossy, 
expensive,  and  heavy  (important  for  air-home  radar  applications). 

•  With  certain  configurations,  less  transmitter  power  is  required  compared  to  the 
monostatic  radar. 

•  Higher  pulse  repetition  frequencies  (PRFs)  can  be  used  because  a  bistatic  system  does 
not  suffer  the  same  range  blindness  as  the  equivalent  monostatie  system. 

•  If  the  target  angle  ean  be  measured  at  both  sites,  as  well  as  the  bistatic  range,  data  ean 
be  cheeked  for  self-consisteney  to  remove  false  alarms. 
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Taking  a  closer  examination  of  the  military  applieations,  bistatie  radar  prineiples 
have  been  deployed  in  systems  such  as  semiactive  homing  missiles,  forward-seatter 
fenees,  multistatie  radar,  and  hitchhikers  [3].  Mueh  research  is  being  done,  and 
eontinues,  in  partieular  with  regard  to  the  last  two  applieations  mentioned. 

The  multistatie  radar  uses  several  distributed  reeeivers  assoeiated  with  a  single  or 
multiple  transmitters.  It  is  suitable  for  wide  area  traeking  and  it  also  inereases  the 
probability  of  deteetion.  The  bistatic  hitchhiker  uses  a  transmitter  of  opportunity,  either 
friendly  (eooperative)  or  hostile  (non-eooperative),  to  deteet  and  loeate  targets  near  the 
transmitter  or  reeeiver  site.  There  is  no  doubt  that  these  bistatie  radar  prineiples  are 
playing  a  bigger  role  in  today’s  military  for  the  deteetion  of  low  signature  targets  and 
improving  the  eleetronie  eountermeasures  and  eounter-eountermeasures. 

Design,  development,  and  testing  of  radar  systems  have  traditionally  been 
expensive  and  time  eonsuming  tasks.  Radar  designers  are  often  faeed  with  many  design 
iterations  to  optimize  the  system  performance.  Using  eomputer  simulation  and  modeling 
provides  the  designer  with  the  essential  tools  to  examine  design  parameters  and  evaluate 
the  system  performanee  before  the  aetual  implementation  in  hardware. 

There  are  many  software  applications  available  for  monostatie  radar  simulation 
but  none  exits  in  the  publie  domain  for  bistatie  radar.  The  objeetive  of  this  thesis  was  to 
develop  a  general  software  model  to  analyze  the  bistatie  radar  system  performanee  as  a 
funetion  of  the  system  design  parameters. 

The  thesis  foeused  on  the  design  and  development  of  a  Graphieal  User  Interfaee 
(GUIs)  in  MATLAB  to  integrate  the  bistatie  radar  parameters,  and  examine  several 
eritieal  performanee  tradeoffs.  The  end  produet  is  a  program  named  Bidarsa,  whieh 
provides  a  eonvenient  instructional  tool  for  engineering  eommunities  to  understand, 
analyze,  and  prediet  bistatie  radar  performanee. 

Chapter  II  provides  the  theoretieal  baekground  required  in  the  formulation  of  the 
software  model.  It  presents  an  overview  of  bistatie  radar  definitions  and  the  parameters 
involved  in  the  range  equation.  The  radar  deteetion  eontour,  the  ealeulation  to  determine 
target  loeation,  eonsiderations  for  the  bistatie  footprint,  the  doppler  relationship  for  a 


2 


moving  target,  and  the  eleetronic  countermeasures  (ECM)  issues  are  discussed  with  their 
supporting  mathematical  equations. 

Chapter  III  presents  the  design  approach  for  the  Bidarsa  program.  A  system 
block  diagram  detailing  the  interfaces  and  data  flow  sequences  is  used  to  illustrate  the 
functions  of  the  eight  major  modules  and  how  they  are  integrated. 

Chapter  IV  uses  a  set  of  radar  system  parameters  compiled  to  demonstrate  the 
Bidarsa  program  functions  and  discuss  its  results. 

Chapter  V  uses  the  Bidarsa  software  to  predict  the  system  performance  in 
different  practical  scenarios  suggested  in  References  [2]  and  [4].  The  result  is  analyzed 
and  compared  with  the  monostatic  radar  configuration. 

Chapter  VI  summarizes  the  research  and  suggests  further  work  that  could  be 
performed  to  enhance  the  code. 

Finally,  the  appendix  lists  the  software  codes  written  for  the  Bidarsa  program. 
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II.  BISTATIC  RADAR  THEORY 


A,  DEFINITIONS 

The  IEEE  Standard  686-1997  [5]  has  defined  bistatic  radar  as  a  radar  using 
antennas  for  transmission  and  reception  at  sufficiently  different  locations  that  the  angles 
or  ranges  from  those  locations  to  the  target  are  significantly  different,  figure  2.1  is  used 
to  illustrate  this  definition,  showing  a  transmitter  Tx  and  receiver  Rx  being  situated  at  two 
locations  with  a  baseline  separation  distance  L.  The  bistatic  angle,  P,  is  one  of  the 
important  parameters  which  characterizes  bistatic  radar  and  affects  system  performance. 
These  relationships  will  be  discussed  in  the  subsequent  sections. 


figure  2.1  Bistatic  radar  geometry. 


There  is  no  strict  definition  specifying  the  distance  between  the  transmitter  and 
receiver  sites.  Several  discussions  [1,  3,  6]  suggest  that  the  baseline,  L,  could  range  from 
tens  of  meters  to  hundreds  of  kilometers,  depending  on  the  required  functions  of  the 


5 


particular  bistatic  radar  system.  The  distanee  Rj  is  the  range  from  the  transmitter-to- 
target,  while  Rr  is  the  range  from  the  target-to-reeeiver. 

B,  RANGE  EQUATION 

The  bistatie  range  equation  is  derived  in  a  similar  manner  as  the  monostatie  radar 
range  equation.  The  bistatie  maximum  range  produet  is  given  by  [3] 

P,G,G,l^c,F}Fl  1’  p 

or 

)„„=*•  (2.2) 

where 

Rt  =  transmitter-to-target  range  (m), 

Rr  =  target-to-reeeiver  range  (m), 

Pt  =  transmitted  power  output  (W), 

Gt  =  transmitting  antenna  power  gain, 

Gr  =  reeeiving  antenna  power  gain, 

X  =  transmitting  frequeney  wavelength  (m), 

Ob  =  target  bistatie  radar  eross  seetion  (RCS)  (m  ), 

Ft  =  pattern  propagation  factor  for  the  transmitter-to-target  path, 

Fr  =  pattern  propagation  faetor  for  the  target-to-reeeiver  path, 

k  =  Boltzmann’s  eonstant,  1.3807  x  10'^^  J/K, 

Ts  =  reeeiver  noise  temperature  (°K), 

Bn  =  reeeiver  noise  bandwidth  (Hz), 

{S/N)rnin  =  minimum  signal-to-noise  power  ratio  required  for  detection. 
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Lt  =  transmitting  system  losses  (>  1), 

Lr  =  reeeiving  system  losses  (>  1), 

K  =  bistatic  maximum  range  product. 

In  the  bistatic  radar  range  equation,  it  is  observed  that  the  maximum  range  {RfRid^ 
has  replaced  the  range  R!*  in  the  monostatic  radar  equation,  due  to  the  different  locations 
of  the  transmitter  and  receiver. 

The  pattern-propagation  factors  for  the  transmitter-to-target  and  target-to-receiver. 
Ft  and  Fr,  also  need  to  be  taken  care  of  separately  in  the  bistatic  case.  This  is  because 
the  transmitting  path  and  receiving  path  can  be  very  different  in  nature  as  compared  to 
the  monostatic  implementation. 


C.  OVALS  OF  CASSINI 


The  oval  of  Cassini  is  the  locus  of  a  points,  the  product  of  whose  distance  from 
two  fixed  points  is  a  constant.  In  the  bistatic  radar  case,  it  is  interesting  to  analyze  the 
system  performance  by  plotting  the  ovals  of  Cassini  as  a  function  of  signal-to-noise 
power  ratio,  S/N,  with  the  ranges,  Rt  and  Rr,  at  the  vertices  of  the  triangle. 

Equation  (2.1)  is  modified  and  written  as 


{RtRr) 


2 

max 


K 


(2.3) 


where  K  is  the  bistatic  radar  constant  and  given  by 

(44kT,B,LrL. 


(2.4) 


The  constant  K  is  related  to  the  bistatic  maximum  range  product,  K ,  by 

K  =  (2.5) 

Figure  2.2  shows  Rt  and  Rr  being  converted  into  the  polar  coordinate  system 
{r,6).  From  geometry 
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Rj^  =  +  1}  !  A^  +  rL  cos  0 


(2.6) 


and 


+1}  ! A^-rLcosG . 


(2.7) 


Combining  Equations  (2.6)  and  (2.7)  yields 

[R^R^f  =(r^  !  A)^  -r^L^  cos'^  0 . 


(2.8) 


Figure  2.2 


Bistatic  radar  geometry  converted  to  the  polar  coordinate  system.  (After 

Ref  [3].) 


The  constant  SIN  ratio  can  subsequently  be  plotted  as  ovals  of  Cassini  using  the 
equation 


(5/A)  = - ^ - . 

(r^+L^IA)  -r^L^cos^G 


(2.9) 


Reference  [3]  has  provided  an  example  of  the  oval  of  Cassini  plot,  as  shown  in  Figure 
2.3,  with  K  arbitrarily  set  to  30F'^  and  constant  SIN  ratios  from  10  dB  to  30  dB. 
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For  analytical  purposes,  erg,  Ft  and  Fr  are  assumed  to  be  invariant  with  r  and  6  so 
as  to  generate  a  constant-SNR  oval-of-Cassini  plot.  A  series  of  plots  shows  that  bistatic 
radar  operation  can  be  divided  into  three  distinct  regions; 

1 .  L>  2-sIlc  .  Two  separate  ellipses  result  enclosing  the  transmitter  and 

receiver.  It  is  considered  receiver  centered  whenR^.  transmitter 

centered  wheni?^  »  Rj. . 

2.  L<  2-s[k  .  A  single  continuous  ellipse  results. 

3.  L  =  2yfK  .  A  lemniscate  with  a  cusp  at  origin  results. 


SVN  •  to  08 


Figure  2.3  Ovals  of  Cassini,  with  constant  SNR  plots  for  K=30L‘^.  (From  Ref.  [3].) 


D.  TARGET  LOCATION 

A  bistatic  receiver  usually  measures  the  angles  of  arrival  of  the  target  echo  and 
estimates  the  range.  The  receiver-to-target  range,  Rr,  eannot  be  measured  direetly,  but 
can  be  calculated  by  solving  the  bistatic  triangle.  When  the  transmitter  location  is  made 
known  to  the  receiver,  the  baseline  range,  L,  is  known.  The  total  signal  propagation  time 
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from  transmitter  to  target  to  receiver  is  also  required  to  be  established.  This  is  needed  for 
the  receiver  to  convert  the  time  measurement  into  the  range  sum  estimation,  Rt  +  Rr. 

From  Figures  2.2  and  2.3,  it  can  be  seen  that  the  isorange  contour  has  to  be  an 
ellipse,  rather  than  the  tangent  approximation.  It  is  defined  by 

Rj,+Rj^=2ci  (2.10) 

where  2a  is  the  major  axis  of  the  ellipse,  the  eccentricity  of  the  ellipse,  e,  is 

e  =  L/2a  (2.11) 

e  =  L/{R^  +  Rj^).  (2.12) 

Using  the  law  of  cosines,  Rr  and  Rt  can  be  calculated  by  measuring  the  L,  (Rt  +  Rr),  and 
either  the  angle  6t  or  Or. 

For  the  case  when  Or  is  measured 

Rl=I^+L^-  2R^L  cos  (90°  -  0j ) 

^  (Rt+rJ-l- 

^  2[Rj  + Rj^-LsmOj) 

1 

Rj^  =  [rI  +  1}  - 2Rj.L sin . 

Combining  Equations  (2.12)  and  (2.14)  yields 

^  2e(l-esin0j,) 

and  combining  Equations  (2.12)  and  (2.15)  yields 

L{^e^  +  \-2R  sin^^) 

^  2e[\-esmOj) 

Similarly,  for  the  case  when  Or  is  measured 

Rl=Rl+L^-  2R^L  cos  (90°  +  ) 


(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 
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2[Rj,  +  Rj^+ LsinOi^) 


(2.19) 


R 


Rj  =  [rI  +I}  +  2Rj^L sin Oj^Y  . 
Combining  Equations  (2.12)  and  (2.19)  yields 

*  2e(l  +  esin^^) 

Finally,  combining  Equations  (2.12)  and  (2.20)  yields 

+  l  +  2i?sin^^) 

^  2e(l  +  esin0^) 


(2.20) 


(2.21) 


(2.22) 


Using  the  law  of  sines,  it  is  also  possible  to  establish  the  parameter  ratios  to 
calculate  the  bistatic  angle,  P,  with  the  measured  Rr  and  Rt  and  the  calculated  Or  and  Or. 
Their  relationship  is  given  by 


Thus 


Also, 


Rj^  Rj  L 

sin(9O°-0j,)  sin  (90°  +  ^^)  sm/3 


(2.23) 


CO?,0j 

(2.24) 

Rj 

cosOg^ 

cosOj 

(2.25) 

L 

sm/3 

sin“' 

^  L  cos  0j  ^ 

(2.26) 

\  J 

Rj. 

cosOg^ 

(2.27) 

L  ~ 

sin  P 
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/5  =  sm  ^ 


(2.28) 


^  L  cos  0^  ^ 

y  Rj  ) 


E.  COVERAGE  AREA 

There  are  two  eonsiderations  for  establishing  the  eoverage  area  of  the  bistatie 
radar:  (1)  the  deteetion-eonstraint  whieh  is  determined  by  the  maximum  range  oval  of 
Cassini,  and  (2)  the  line-of-sight  (LOS)  eonstraint  whieh  is  determined  by  the  geometry 
between  the  transmitter,  target,  and  reeeiver. 


1.  Detection-Constrained  Coverage 

The  eoverage  is  determined  by  the  area  within  a  maximum-range  oval  of  Cassini, 
whieh  is  established  by  the  bistatie  maximum  range  produet,  /r,  in  a  bistatie  radar 
system. 

Depending  on  region  of  operation,  i.e.,  the  relationship  between  L  and  24k:  , 
different  expressions  for  the  area,  Ab,  apply.  The  detailed  derivation  of  eaeh  expression 
is  provided  in  Appendix  D  of  [3]. 

a.  L<  24k:  .  A  single  eontinuous  ellipse. 


A,,=7rK 


1- 


\6k^ 


1 

viy 


1»3 

V2.4y 


r 


\6k^ 


rn 

UJ 


(2.29) 


Ag^  ~  7CK 


1- 

f  1  ^ 

164, 

y  ) 

U6384j 

(2.30) 


b.  L>  244  .  Two  separate  ellipses  enelosing  the  transmitter  and  reeeiver. 


^B1 


27CK^ 


1} 


1+ 


-I- 


f  1  ^ 

2 

f  16r"^ 

-1- 

f  3"  ^ 

2 

( \6k^'" 

1  ^  ; 

[2".3!*2!  y 

\  ^  J 

2 

[2^41.3!  y 

\  ^  J 

-!-• 


(2.31) 
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(2.32) 


A 


iTtK^ 


B2 


\  2k^  12r^  lOOr" 

1  +  — 


V 


n 


n 
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c.  L  =  I'Jk  .  A  lemniscate  with  a  cusp  at  the  origin 

A^,=2k.  (2.33) 

Comparing  a  bistatic  radar  with  maximum  range  produet,  /r,  to  a  monostatie 
radar  with  an  equivalent  monostatie  range,  Vr  ,  it  is  found  that  the  area  of  eoverage  for  a 
bistatie  radar  is  smaller  than  that  of  the  monostatie  radar  [3]. 


2.  Line-of-Sight-Constrained  Coverage 

The  bistatie  line-of-sight  (LOS)  eoverage  is  affeeted  by  multipath,  diffraetion, 
refraction,  and  shadowing,  including  the  earth  curvature.  The  coverage  area,  Ac,  is 
established  by  the  eommon  area  eovered  by  both  the  transmitter  and  receiver  eoverage 
eireles,  as  shown  in  Figure  2.4.  The  4/3  earth  model  gives  radii  of  eoverage  eircles  as 
follows 

=130(7^  +  7^)  (2.34) 


(2.35) 


where 

ht  =  target  altitude  (km), 

hr  =  transmitter  antenna  altitude  (km), 

Hr  =  reeeiver  antenna  altitude  (km). 

For  L  +  Tj,  >  rj^>  L-Tj.  or  L  +  rj^>  Tj,  >  L-Tj^,  the  interseetion  of  the  two  eoverage 
eireles  determines  the  eommon  coverage  area 

=  I  [<I>R  -  sin  (l>R  )  +  r^[(l)T-  sin  (2.36) 
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where 


^J=l  cos 


-1 


A  2  1  , 

Vj.  +L 

ItjL 


A  2  2  ,  r2  A 

,  -  _1  -Ty.  +L 

=  2  cos  ^ 


V 


2r^Z 


y 


otherwise, 


if  rj.  +  rj^<L,  then  =  0; 


(2.37) 

(2.38) 


(2.39) 


if  rj.>L  +  rjj,  then  ; 


(2.40) 


if  ,  then  =  TTr^ . 


(2.41) 


Figure  2.4  Coverage  area  geometry.  (From  Ref  [3].) 

F.  CUTTER  CELL  AREA 

The  intersection  of  the  range  cell  and  the  mainbeam  footprint  is  used  to  determine 
the  bistatic  mainlobe  clutter  cell  area,  Ac.  It  is  useful  to  determine  if  a  range  cell  is 
beamwidth  limited  or  pulsewidth  limited. 
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Figure  2.5  illustrates  the  beamwidth  limited  case  with  clutter  cell  area  given  by 


[Rj^A9r)[RtA0t)  ^2  42) 

sm/3 

The  terms  andi?^A^^  are  the  cross-range  dimension  of  the  transmitting  and 

receiving  beam  at  the  clutter  cell,  respectively.  The  angles  AO^  and  A9r  are  the  3-dB 
beamwidths  of  the  transmitting  and  receiving  antennas,  respectively. 


Figure  2.5  Beamwidth  limited  clutter  cell  area.  (From  Ref.  [7].) 


Figure  2.6  illustrates  the  pulsewidth  limited  case  with  clutter  cell  area  given  by 


COS 


(A/2) 


ctRrAOr 

2cos^  {^/2) 


(2.43) 


The  term  t  is  the  compressed  pulsewidth,  whereas  AR  is  the  approximation  to  the 
bistatic  range  cell  given  by 
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CT 


(2.44) 


2cos[P  H) 


TX 


Figure  2.5  Pulsewidth  limited  clutter  cell  area.  (From  Ref.  [7].) 


Appendix  B  of  Reference  [3]  has  provided  an  estimation  for  the  maximum  error 
of  a  range  cell,  .  The  maximum  error  occurs  on  the  perpendicular  bisector  of  the 
baseline,  where  is  a  maximum.  It  is  given  by 


a{a  -a) 
~  b{b'-b) 


(2.45) 


where 


a  =  semimajor  axis  of  inner  ellipse  =  {Rj  +  7?^  )  /  2 , 
a  =  semimajor  axis  of  outer  ellipse  =  (a  +  cf/  2) , 

b  =  semiminor  axis  of  inner  ellipse  =  {a^  -l}  !  Ay  , 

(af -ZV4? 


b'  =  semiminor  axis  of  outer  ellipse  = 
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G.  MAXIMUM  UNAMBIGUOUS  RANGE  AND  PRF 

The  maximum  unambiguous  range  for  a  bistatic  radar  is  given  by 


[Rj,  +  Rj^)^ 


c 

PRF 


and  the  maximum  unambiguous  PRF  is  given  by 


c 

{Rj.+Rj^)^^ 


c 

+  2/r(l  +  coSy5) 


(2.46) 


(2.47) 

(2.48) 


H.  DOPPLER  RELATIONSHIP 


The  doppler  shift,  fs,  is  the  time  rate  of  change  of  the  total  path  length  of  the 
scattered  signal,  normalized  by  the  wavelength.  For  bistatic  radar,  the  doppler  shift  is 
given  by 


(2.49) 


/.= 


1 

I 


dRr 

— -+ 

dt 


dRg^ 

dt 


(2.50) 


Figure  2.6  shows  a  bistatic  doppler  geometry.  The  target’s  velocity  vector  has  magnitude 
V  and  aspect  angle  d.  They  are  being  projected  to  the  transmitter  and  receiver  as  Vj,  Vr 
and  dr,  dR  respectively.  A  stationary  transmitter  and  receiver  will  have  the  projection  of 
the  target  velocity  vector  onto  the  transmitter-to-target  LOS  and  receiver-to-target  LOS 
as  follows 


^  =  Vcos{d-/5l2)  (2.51) 


-^  =  Vcos{d  +  pl2).  (2.52) 

dt 
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Combining  Equations  (2.50),  (2.51)  and  (2.52)  yields 

f^  =  {V !  X)[co^{d-l3l2)  +  co^[d+pl2)]  (2.53) 

fs=[2V !  X)cos5cos[l3l2).  (2.54) 


I.  ELECTRONIC  COUNTERMEASURES 


Bistatic  radar  can  also  be  subjected  to  intentional  interference  initiated  by  hostile 
transmissions  from  other  sources.  Thus  it  is  essential  to  examine  how  the  performance 
will  be  affected  under  such  an  environment.  When  a  noise  jammer  is  employed,  the  radar 
receiver  will  be  subjected  to  the  jamming  power  spectrum  density,  Jo,  given  by 


(2.55) 


where 


Pj  =  jammer  transmitted  power  (W), 
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Gj  =  j  ammer  transmit  antenna  power  gain, 

Gr  =  radar  reeeiving  antenna  power  gain, 

Fj  =  j  ammer  pattern  propagation  faetor  =  Fjfj  , 

Fj  =  propagation  faetor, 

fj  =  jammer  antenna  patter  faetor,  refereneed  to  the  reeeiver  site, 

/rj  =  radar  reeeiving  antenna  patter  faetor,  refereneed  to  jammer  site, 

Bj  =  j  ammer  bandwidth  (Hz), 

Rj  =  j  ammer- to-reeeiver  range  (m), 

Lj  =  jammer  system  losses  (>  1). 

Assuming  the  ease  of  a  single  jammer,  substituting  Equation  (2.55)  into  (2.1),  the  bistatie 
bum-through  equation  is  given  by 


{RtRr  ) ,, 


PG  T  B 


P.,G,, 


F 


LtLr 


Bn  47r(5/A) 


(2.56) 


J.  SUMMARY 

The  bistatie  radar  equation  is  derived  in  a  manner  similar  to  the  monostatie  radar 
range  equation,  with  the  bistatie  range  (RtFrY  replaeing  the  range  R'^  in  the  monostatie 
radar  ease.  A  eonstant  SNR  plot  has  the  shape  of  an  oval  of  Cassini.  By  resolving  the 
bistatie  triangle,  it  is  possible  to  ealeulate  the  target  loeation.  The  radar  eoverage  is 
defined  by  the  region  where  the  target  is  deteetable  by  the  reeeiver  and  within  the  LOS  of 
both  transmitter  and  reeeiver.  The  bistatie  footprint  ean  be  either  beamwidth  or 
pulsewidth  limited.  Both  doppler  relationships  and  eleetronie  eountermeasures 
eonsiderations  were  also  diseussed.  Chapter  III  uses  the  theory  diseussed  in  this  ehapter 
to  establish  the  approaeh  for  the  bistatie  radar  simulation,  and  the  integration  of  the 
various  software  modules. 
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III.  SOFTWARE  DESIGN  APPROACH 


A,  OVERVIEW 

Bidarsa  is  a  bistatic  radar  simulation  developed  using  the  graphieal  user  interfaees 
(GUIs)  available  in  MATLAB.  The  program  is  designed  to  provide  easy  user  interfaees 
for  radar  designers,  engineers,  and  students  to  examine  the  various  design  parameters  for 
a  bistatie  radar  and  analyze  its  performanee. 

The  bistatie  radar  model  bloek  diagram  is  shown  in  Figure  3.1.  Eight  main 
modules  were  developed  to  earry  out  different  bistatie  radar  system  analyses  as  diseussed 
in  Chapter  II.  Figure  3.1  also  indieates  the  data  flow  sequenees.  It  illustrates  how  eertain 
modules  will  need  to  be  exeeuted  prior  to  the  others,  as  their  results  are  used  in 
subsequent  analyses  by  the  other  modules. 

The  Bidarsa.m  module  allows  a  user  to  seleet  one  of  the  three  group  funetions:  (1) 
speeify  the  system  parameters,  (2)  earry  out  a  basie  system  performanee  analysis,  or  (3) 
to  earry  out  an  extended  system  performanee  analysis. 

The  Defpara.m  module  provides  a  eonvenient  way  for  a  user  to  speeify  the 
bistatie  radar’s  transmitter  and  reeeiver  design  parameters.  It  also  allows  target 
information  and  the  environmental  information  to  be  speeified.  The  parameters  ean  be 
updated  at  any  point  in  time  to  allow  a  new  seenario  to  be  examined,  and  its  effeet  on  the 
overall  system  performanee  evaluated. 

The  RngCal.m  module  eomputes  the  bistatie  maximum  range  produet,  bistatie 
radar  eonstant,  and  the  maximum  range-unambiguous  pulse  repetition  frequeney  (PRF). 
This  module  allows  the  user  to  earry  out  system  analysis  on  four  different  design 
parameters:  (1)  transmitter  power,  (2)  system  frequeney,  (3)  target  RCS,  and  (4)  signal- 
to-noise  ratio.  The  user  is  also  able  to  plot  the  deteetion  eontour  and  examine  the  oval  of 
Cassini  with  four  different  eonstant  SNR  plots. 
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Figure  3.1  Bistatic  radar  system  analysis  model  block  diagram. 
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The  TgtLoc.m  module  uses  the  range-sum  value  and  angle  6  to  eompute  the 
loeation  of  the  target.  It  provides  the  transmitter-to-target  and  target-to-reeeiver  ranges 
and  ealeulates  the  bistatie  angle. 

The  CovAre.m  module  evaluates  the  eoverage  areas  for  both  deteetion-eonstrained 
and  LOS-eonstrained  oases.  The  module  also  provides  an  analysis  on  the  transmitter- 
target-reoeiver  height  relationship  requirement,  and  the  effeots  of  different  baseline 
ranges. 

The  FotPrt.m  module  uses  the  transmitter- target-reoeiver  ranges,  antenna 
beamwidth  and  transmitter  pulsewidth,  to  evaluate  the  bistatie  footprint  and  olutter  area. 
The  module  examines  the  results  and  determines  if  the  deteotion  is  beamwidth  or 
pulsewidth  limited.  The  range  oell  size  and  maximum  range  errors  are  also  oomputed  in 
this  module. 

The  DopRel.m  module  examines  the  doppler  shift  effeots  of  a  moving  target.  It 
also  provides  an  analysis  of  the  target-velooity  aspeot  angle  versus  bistatie  angle,  and 
how  they  affeot  the  doppler  frequenoy  shift. 

The  EWEff.m  module  provides  a  means  to  evaluate  the  bistatie  radar  performanoe 
under  different  jammer  environments.  The  module  provides  an  analysis  on  the  bum- 
through  range  for  different  target  RCS  and  jammer  power  soenarios. 

A  Help  menu  is  provided  in  all  modules  to  assist  the  user  to  navigate  and 
understand  how  to  use  eaoh  module. 

Print  and  Close  funotions  are  also  designed  into  all  modules  to  allow  the  user  to 
make  a  hard  eopy  and  exit  the  seleeted  menu  whenever  desired. 

The  detailed  MATLAB  programs  for  the  eight  modules  are  listed  in  the  Appendix 
of  this  thesis. 
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B. 


BIDARSA  MAIN  MENU 


To  run  the  Bidarsa  GUI,  the  user  has  to  start  MATLAB  and  set  the  path  of 
direetory  to  where  the  Bidarsa  eode  is  stored.  Typing  Bidarsa  at  the  eommand  window 
will  exeeute  the  GUI  as  shown  in  Figure  3.2.  Among  the  three  groupings,  the  user  should 
define  the  system  parameters  first  before  running  the  other  analysis  funetions.  An 
explanation  of  the  Bidarsa  GUI  may  be  obtained  by  elieking  on  the  Help  button. 


Figure  3.2  main  menu  GUI. 
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c. 


DESIGN  PARAMETERS  MENU 


This  module  is  designed  to  collect  the  system  parameters,  target  information  and 
environmental  data.  These  inputs  are  generally  required  by  all  the  other  modules.  It 
functions  as  a  central  database  for  the  program,  and  converts  each  parameter  to  a  suitable 
unit  for  subsequent  computations.  A  pop  up  window  is  designed  to  inform  the  user  that 
new  parameters  are  being  loaded  for  any  follow-on  analysis.  Figure  3.3  shows  the  GUI 
layout  designed  for  the  module.  A  typical  value  for  each  parameter  was  chosen  as  the 
default  value  for  demonstration  and  initial  calculation  purposes. 


Figure  3.3  Define  parameter  GUI. 
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D.  RANGE  CALCULATION  MENU 


This  module  is  designed  to  earry  out  the  most  fundamental  analysis  for  a  bistatic 
radar  system.  Several  other  modules  use  its  result  to  carry  out  further  analysis.  Figure 
3.4  shows  the  GUI  layout  designed  for  the  module.  With  the  input  data  from  the 
Defpara.m  module,  it  computes  the  following: 

1 .  Bistatic  maximum  range  product,  {RjRn  ,  as  detailed  in  Equation  (2.1). 

2.  Bistatic  radar  constant,  K,  as  detailed  in  Equation  (2.4). 

3.  Maximum  range-unambiguous  PRE,  as  detailed  in  Equation  (2.48). 

4.  System  Analysis  for 

a.  Transmitting  power  selection 

b.  System  frequency  selection 

c.  Target  RCS  variation 

d.  Receiver  SNR  variation 

5.  Detection  contour  plot  for  ovals  of  Cassini,  as  detailed  in  Equation  (2.9). 


Eigure  3.4  Range  calculation  GUI. 
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E. 


TARGET  LOCATION  MENU 


This  module  is  designed  to  estimate  the  location  of  the  target  using  the  range-sum 
value  and  measured  angle  9.  Figure  3.5  shows  the  GUI  layout  designed  for  the  module. 
Only  one  of  the  angles,  9t  or  Or,  is  required  for  the  reference,  as  discussed  in  Chapter  II 
Section  D.  It  computes: 


1 .  Transmitter-to-target  range,  as  detailed  in  Equations  (2.16)  and  (2.22). 

2.  Target-to-receiver  range,  as  detailed  in  Equations  (2.17)  and  (2.21). 

3.  Bistatic  beta  angle,  as  detailed  in  Equations  (2.26)  and  (2.28). 


Eigure3.5  Target  location  GUI. 
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F.  COVERAGE  AREA  MENU 

This  module  is  designed  to  evaluate  the  region  or  area  on  the  bistatie  plane  where 
the  target  is  detectable  by  the  receiver  and  within  the  LOS  of  both  transmitter  and 
receiver.  The  user  can  specify  a  different  altitude  for  the  target,  transmitter  antenna,  and 
receiver  antenna  for  this  analysis.  Figure  3.6  shows  the  GUI  layout  designed  for  the 
module.  It  computes: 

1.  The  region  of  operation,  based  on  the  relationship  of  L  and  2-Jk ,  as 
discussed  in  Chapter  II  Section  E. 

2.  Coverage  area  for  the  detection-constrained  case,  as  detailed  in  Equations 
(2.30),  (2.32),  and  (2.33). 

3.  Coverage  area  for  the  EOS-constrained  case,  as  detailed  in  Equations 
(2.36),  (2.39),  (2.40),  and  (2.41). 

4.  System  analysis  for 

a.  Transmitter  antenna  altitude  requirements. 

b.  Receiver  antenna  altitude  requirements. 

c .  T arget  altitude  requirements . 

d.  Baseline  distance  variations. 
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Figure  3.6  Coverage  area  GUI. 
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G. 


BISTATIC  FOOTPRINT  AND  CLUTTER  AREA  MENU 


This  module  is  designed  to  evaluate  the  bistatic  footprint  and  clutter  area  of  the 
bistatic  radar  system.  The  user  can  specify  the  transmitter-target-receiver  ranges,  antenna 
beamwidth  and  transmitter  pulse  information  for  different  analysis.  Figure  3.7  shows  the 
GUI  layout  designed  for  the  module.  It  computes: 

1 .  Whether  the  system  is  a  beamwidth  or  pulsewidth  limited  case. 

2.  The  bistatic  footprint  area,  as  detailed  in  Equations  (2.42)  and  (2.43). 

3.  The  range  cell  size,  as  detailed  in  Equation  (2.44). 

4.  Maximum  error  in  the  range  cell,  as  detailed  in  Equation  (2.45). 


F igure  3 . 7  Clutter  area  GUI . 
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H,  DOPPLER  RELATIONSHIP  MENU 


This  module  is  designed  to  evaluate  the  doppler  shift  effects  of  a  moving  target. 
The  user  can  provide  different  target  velocities  and  aspect  angles  for  the  analysis.  Figure 
3.8  shows  the  GUI  layout  designed  for  the  module.  It  computes: 

1 .  Bistatic  doppler  shift,  as  detailed  in  Equation  (2.54). 

2.  System  analysis  for 

a.  Target  velocity  aspect  angle  variation. 

b.  Bistatic  angle  variation. 

c.  Bistatic  doppler  shift  variation. 


Figure  3.8  Doppler  relation  GUI. 
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I.  EW  EFFECTS  MENU 

This  module  is  designed  to  evaluate  the  bistatic  radar’s  performance  when 
subjected  to  a  jamming  transmission.  The  user  can  define  the  jammer  specifications  for 
simulating  different  operational  scenarios.  Figure  3.9  shows  the  GUI  layout  designed  for 
the  module.  It  computes: 

1.  Bistatic  maximum  range  product  with  jammer,  as  detailed  in  Equation 

(2.56). 

2.  System  analysis  for 

a.  Jammer  power  variation. 

b.  Target  RCS  variation. 

c.  Burn-through  range  analysis. 


Figure  3.9  EW  Effects  GUI. 
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J. 


SUMMARY 


The  MATLAB  GUI  is  used  to  develop  the  Bidarsa  simulation  program,  as  it 
offers  easy  user  interfaees  and  excellent  data  visualization.  Eight  modules  were 
developed  and  integrated  as  a  complete  suite.  The  Bidarsa  program  provides  a  user  the 
options  to  specify  the  bistatic  radar  system  parameters,  the  target  information  and 
environmental  data.  Several  critical  system  performance  measures  are  calculated,  and 
data  presentation  tailored  to  the  study  of  design  tradeoffs.  The  bistatic  doppler 
relationship  and  electronic  countermeasures  considerations  are  also  included  to  provide 
extended  analyzes.  Chapter  III  provides  sample  analyzes  on  a  set  of  radar  parameters. 
Plots  generated  from  the  different  modules  will  also  be  discussed. 
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IV.  DESCRIPTION  OF  SOFTWARE  OPERATION 


A,  TEST  DATA 

To  evaluate  the  Bidarsa  program,  a  set  of  system  parameters  were  eompiled  from 
various  published  sourees  [1,3]  and  used  for  a  detailed  system  analysis.  The  test  values 
are  as  shown  in  Table  1. 


Table 


Transmitting  power 

50dBW 

Transmitter  antenna  gain 

30  dB 

Frequeney 

10  GHz 

Reeeiver  min  SNR 

0  dBm 

Reeeiver  antenna  gain 

30  dB 

Reeeiver  noise  temperature 

300  K 

Reeeiver  noise  bandwidth 

60  MHz 

Baseline  range 

100  km 

Target  RCS 

10  dBm^ 

Bistatie  beta  angle 

1  degree 

Transmitter  system  loss 

1  dB 

Receiver  system  loss 

1  dB 

Tx-Tgt  pattern  propagation  factor 

-3.5  dB 

Tgt-Rx  pattern  propagation  factor 

-3.5  dB 

Test  parameters  used  to  evaluate  Bidarsa  program. 


35 


The  Bidarsa  program  is  run  and  the  above  values  are  entered  and  updated  in  the  Defpara 
menu.  The  results  of  all  module  caleulations  are  displayed  in  the  MATLAB  Command 
Window. 

B,  BISTATIC  RANGE  PRODUCT  ANALYSIS 

Selecting  the  Range  Calculation  menu,  the  program  computed  the  following 
results  for  the  three  range  related  parameters,  as  shown  in  Figure  4.1.  Selecting  one  of 
the  four  system  analysis  buttons  on  the  menu  will  generate  the  respective  plot  as  shown 
in  Figure  4.2.  When  the  plot  detection  contour  button  is  selected,  the  program  generates 
a  plot  for  the  oval  of  Cassini  for  constant  SNR  at  93  dB,  96  dB,  98.65  dB  and  102  dB,  as 
shown  in  Figure  4.3. 


Bidarsa  -  Range  Calculation: 

Bistatic  Maximum  Range  Product  RtRr  is  6.7706e+006  km^. 

Bistatic  Radar  constant  K  is  4.5841e+016. 

Maximum  Range-Unambiguous  PRF  is  1.8299e+000  kHz. 

Figure  4. 1  Range  calculation  menu  output. 

From  Equations  (2.48)  and  (2.46),  it  is  interesting  to  note  that  in  a  bistatic  radar 
system,  increasing  the  bistatic  angle,  p,  will  increase  the  maximum  range-unambiguous 
PRF,  which  will  in  turn,  reduce  the  bistatic  radar  maximum  unambiguous  range. 

Figure  4.3  resembles  Figure  2.3  in  Reference  [3].  With  the  system  parameters 
used  in  this  analysis,  the  lemniscate  {L  =  2'Jk)  occurs  when  the  SNR  is  98.65  dB;  two 
separate  ellipses  enclosing  the  transmitter  and  receiver  {L>  2 Vr  )  are  observed  at  a  SNR 

of  102  dB;  two  single  continuous  ellipses  {L<  2Vr)  are  observed  at  SNRs  of  93  dB  and 
96  dB. 
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Figure  4.2  System  analysis;  transmit  power  (top  left),  frequeney  (top  right),  target 
RCS  (bottom  left)  and  SNR  (bottom  right). 


Reference  [3]  also  suggests  a  few  possible  military  applications  for  each  of  the 
operating  regions. 

1.  For  Z  >  2^/^ ,  with  a  small  oval  around  the  receiver,  i.e.,  the  receiver- 
centered  region,  the  system  can  be  used  for 

a.  Air-to-ground  operations  using  a  stand-off  transmitter  and 
penetrating  aircraft  with  silent  receiver. 

b.  A  semiactive  homing  missile. 


37 


2.  For  L  >  2yjlc ,  with  a  small  oval  around  the  transmitter,  i.e.,  the 
transmitter-centered  region,  the  system  can  be  used  for 

a.  Monitoring  of  activity  near  a  non-cooperative  transmitter. 

3.  For  L  <  2\[k  ,  with  a  single  continuous  ellipse,  the  system  can  be  used  for 

a.  Medium-range  air  defense. 

b.  Missile  tracking  from  ground-based  sites. 


Figure  4.3  Constant  SNR  oval  of  Cassini  plot. 
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C.  TARGET  LOCATION  ANALYSIS 

In  the  Target  Location  menu,  test  values  for  the  range-sum  and  angle  Or  as  shown  in 
Table  2,  are  used  for  the  evaluation.  The  program  eomputed  the  results  for  the  three 
target  loeation  related  parameters,  as  shown  in  Figure  4.4. 


Range-sum  -I-  ) 

350  km 

Angle  Or 

30  degrees 

Table  2.  Test  values  used  for  target  loeation  analysis. 


Bidarsa  -  Target  Loeation: 

Theta  R  Seleeted 

Distanee  of  Target  to  Transmitter 

is  2.0938e+002  km. 

Distanee  of  Target  to  Reeeiver 

is  1.4063e+002  km. 

The  Beta  Angle  ealeulated 

is  24.43  deg. 

Figure  4.4 

Target  loeation  menu  output. 

The  user  ean  toggle  the  angle  6  entry  and  verify  from  the  results  that  the  isorange 
eontour  has  to  be  an  ellipse,  rather  than  the  tangent  approximation.  The  possibility  of 
using  either  Or  or  Or  for  eomputation  offers  more  flexibility  in  the  analysis. 
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D, 


COVERAGE  AREA  ANALYSIS 


In  the  Coverage  Area  menu,  the  test  values  for  the  transmitter,  receiver,  and  target 
altitude  are  shown  in  Table  3.  The  program  computed  the  following  results  in 
determining  the  region  of  operation  and  the  coverage  areas,  as  shown  in  Figure  4.5.  A 
system  analysis  of  the  transmitter-receiver-target  height  requirement  versus  baseline 
range  is  also  generated  as  shown  in  Figure  4.6. 


Table  3. 


Target  height 

30  km 

Transmitter  antenna  height 

10  km 

Receiver  antenna  height 

10  km 

Test  values  used  for  coverage  area  ana^ 


ysis. 


Bidarsa  -  Coverage  Area: 

Target  is  in  the  Single  Cosite  Region. 

2 

The  Detection-constrained  coverage  area  is  2.1271e+007  km  . 

Rx  and  Tx  is  in  LOS. 

The  LOS-constrained  coverage  area 

is  3.7384e+006  km^. 

Figure  4.5  Coverage  area  menu  output. 


In  this  analysis.  Figure  4.6  shows  that,  for  a  baseline  range  T  of  100  km,  an 
adequate  LOS  is  established  for  a  target  height  >  700  m,  only  if  the  transmitting  or 
receiving  antenna  is  >  100  m.  If  the  baseline  is  decreased  to  50  km  (0.5L),  an  adequate 
LOS  is  established  for  a  target  height  >  700  m,  when  the  antenna  is  near  zero  altitude. 
However,  if  the  range  is  increased  to  1200  km  (L2L),  to  achieve  an  adequate  LOS  for 
target  height  >  700  m,  the  antenna  will  need  to  be  raised  to  at  least  180  m. 
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F igure  4.6  T ransmitter-target-rec  eiver  height  analysis . 
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E, 


BISTATIC  FOOTPRINT  ANALYSIS 


In  the  Bistatic  Footprint  &  Clutter  Area  menu,  the  radar  parameters  as  shown  in  Table  4 
are  used  for  the  evaluation.  The  program  determines  if  the  bistatic  footprint  is  limited  by 
beamwidth  or  pulsewidth,  as  shown  in  Figure  4.7.  The  range  cell  size  and  maximum 
error  were  also  generated. 


Receiver-target  range 

100  km 

Transmitter-to-target  range 

50  km 

Receiver  antenna  beamwidth 

10  degrees 

Transmitter  antenna  beamwidth 

10  degrees 

Transmitter  pulsewidth 

30  psec 

Table  4.  Test  parameters  for  bistatic  footprint  analysis. 


Figure  4.7  Bistatic  footprint  menu  output. 


It  can  be  gathered  from  the  Equations  (2.42)  to  (2.45)  that  the  radar  pulsewidth 
plays  a  bigger  role  than  the  beamwidth  in  determining  whether  the  bistatic  footprint  is 
limited  by  beamwidth  or  pulsewidth.  This  effect  can  be  verified  using  this  module. 
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It  is  also  verified  in  this  module  that  variation  in  the  antenna  beamwidth  will  only 
affect  the  beamwidth  limited  bistatic  footprint  area;  it  has  a  small  effect  on  the  range  cell 
size  and  does  not  affect  the  maximum  error  calculated. 

F.  DOPPLER  RELATION  ANALYSIS 

In  the  Doppler  Relation  menu,  the  test  values  shown  in  Table  5  are  used  for  the 
evaluation.  The  program  computes  the  bistatic  doppler  shift  as  shown  in  Figure  4.8.  A 
system  analysis  for  doppler  frequency  shift  versus  different  target  aspect  angle  and 
bistatic  angle  is  generated  as  shown  in  Figure  4.9. 


Target  velocity 

300  m/s 

Target  velocity  aspect  angle 

80  degrees 

Table  5.  Test  parameters  used  for  doppler  relationship  analysis. 


Bidarsa  -  Doppler  Relation: 

Bistatic  Doppler  Shift  is  3.4597e+001  kHz. 

Figure  4.8  Doppler  Relation  menu  output. 

An  analysis  of  Figure  4.9  illustrates  the  following  observations: 

1.  For  all  bistatic  angles,  P,  if  the  target  velocity  aspect  angle,  6,  is  between 
-90°  and  +90°,  there  will  be  a  positive  bistatic  target  doppler. 

2.  For  all  bistatic  angles,  p,  if  the  target  aspect  angle  is  normal  to  the  bistatic 
bisector,  6  =  ±90°,  there  will  be  zero  bistatic  target  doppler. 

3.  For  all  bistatic  angles,  P  <  120°,  if  the  target  aspect  angle  is  collinear  with 
the  bistatic  bisector,  a  maximum  bistatic  target  doppler  will  occur  when 
^  =  0°  and  a  minimum  will  occur  when  d  =  180°  . 
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2jFigure  No.  1 
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Figure  4.9  Doppler  frequency  shift  analysis  with  different  target  aspect  angle  and 

bistatic  angle. 
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G.  ELECTRONIC  COUTNERMEASURES  ANALYSIS 

In  the  EW Ejfects  menu,  the  jammer  specifications  shown  in  Table  6  are  used  for 
the  evaluation.  The  program  computes  the  bistatic  maximum  range  product  with 
jamming  as  shown  in  Figure  4.10.  A  system  analysis  for  the  burn-through  ranges  with 
different  jammer  transmit  powers  and  target  RCS  is  generated  in  Figure  4.11. 


Jammer  transmitted  power 

30dBW 

Jammer  antenna  gain 

5  dB 

Jammer  propagation  factor 

1  dB 

Jammer  bandwidth 

1000  MHz 

Jammer-receiver  range 

100  km 

Jammer  system  loss 

3  dB 

Table  6.  Test  parameters  used  for  electronic  countermeasures  analysis 


Bidarsa  -  EW  Effects: 

Bistatic  Max  Range  Product  with  Jamming  (RtRr)j  is  3.6418e+008  km^. 

Eigure  4.10  EW  Effects  menu  output. 


45 


j^Jpigure  No.  1 
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Burn-through  Range  Analysis 


Figure  4.11  Burn- through  range  analysis  with  different  target  RCS. 


H,  SUMMARY 

A  set  of  radar  parameters  was  used  to  demonstrate  the  simulation  and  analysis 
capability  provided  by  the  Bidarsa  program.  The  simulation  results  from  the  different 
modules  were  discussed  along  with  their  respective  plots.  The  Bidarsa  program  offers 
flexibility  for  the  user  to  study  how  the  individual  parameters  will  affect  the  overall 
system  performance.  Chapter  V  presents  more  examples  to  illustrate  how  the  Bidarsa 
program  can  be  used  to  study  specific  bistatic  configurations  and  scenarios. 
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V.  SAMPLE  SYSTEMS  ANALYSIS 


A,  OVERVIEW 

Chapter  IV  has  shown  the  various  components  of  the  software  simulation  and 
analysis  offered  by  the  Bidarsa  program.  In  this  chapter,  the  use  of  the  Bidarsa  program 
is  further  illustrated  with  two  specific  bistatic  radar  applications. 

B,  MONITORING  OF  SHIPPING  CHANNEL 

Chapter  13  of  Reference  [3]  has  provided  a  scenario  whereby  an  aircraft  uses  a 
bistatic  radar  system  to  monitor  the  shipping  across  a  channel.  A  1-GHz,  1-W 
transmitting  buoy  is  dropped  in  a  shipping  channel  such  that  ships  must  pass  within  1  km 
of  it.  The  aircraft  is  flying  directly  overhead  listening  for  the  echoes  with  a  receiving 
antenna  gain  of  20  dB.  The  buoy  has  a  unity  antenna  gain  and  the  aircraft  receiver 
requires  a  detection  signal  strength  of  -130  dBW.  The  bistatic  RCS  of  a  ship  is  30 
dBsm,  and  all  the  losses  and  multipath  effects  are  neglected. 

The  above  parameters  were  entered  into  the  Bidarsa  program  for  the  analysis. 
The  results  are  tabulated  as  shown  in  Table  7. 


Bistatic  maximum  range  product,  RtRy 

6.7359  km^ 

Bistatic  radar  constant,  K 

4.5372e-006 

Distance  of  target-to-transmitter 

1  km 

Distance  of  target-to-receiver 

6.7359  km 

Beta  angle 

82.17  degrees 

Table  7.  Results  for  shipping-channel 

.  example. 

The  above  bistatic  radar  result  was  used  to  compare  its  performance  to  an  aircraft 
using  monostatic  radar  in  the  same  scenario.  The  transmitting  antenna  is  assumed  to 
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have  a  unity  gain  for  consistency.  The  monostatic  radar  will  give  a  maximum  detection 
range  of  only  2.595  km,  due  to  the  factor  of  in  the  monostatic  radar  range  equation. 
The  bistatic  radar  system  in  this  application  has  a  range  advantage  over  the  monostatic 
radar. 

This  comparison  has  implied  that  by  using  the  bistatic  radar  configuration,  the 
aircraft  could  fly  at  a  much  higher  altitude  to  monitor  the  shipping  channel  without  being 
detected.  On  the  other  hand,  if  the  monostatic  radar  uses  a  higher  gain  antenna,  it  will 
make  up  for  the  monostatic  factor  and  require  less  transmitting  power  for  the  same 
range  performance.  However,  the  monostatic  radar  configuration,  being  an  active 
system,  might  possibly  give  away  the  position  of  the  aircraft. 

The  transmitting  power,  system  frequency,  target  RCS,  and  SNRmin  versus 
maximum  range  product  are  shown  in  Figure  5.1. 


Figure  5.1  System  analysis  for  shipping-channel  example:  transmit  power  (top  left), 
frequency  (top  right),  target  RCS  (bottom  left)  and  SNR  (bottom  right). 
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Depending  on  the  requirements  of  the  aircraft  operation  profile,  Figure  5.1 
illustrates  the  tradeoffs  in  the  system  design  parameters,  so  as  to  achieve  a  higher 
probability  of  detection  and  aircraft  survivability. 

C.  NFS  BISTATIC  RADAR  EXPERIMENT 

An  experiment  was  conducted  at  the  Naval  Postgraduate  School  (NPS)  in  1997  to 
investigate  the  feasibility  of  using  direct  television  broadcast  satellites  and  ground-based 
bistatic  receivers  to  detect  targets  [4].  With  the  Hughes  Direct  Broadcast  Satellite  (DBS- 
1)  used  to  illuminate  the  target,  receive  signals  from  the  direct  path  antenna  and  bistatic 
receive  antenna  were  compared  and  processed  to  detect  the  presence  of  targets.  The 
system  configuration  is  shown  in  Figure  5.2. 

The  system  parameters  used  for  the  experiment  are  as  shown  in  Table  8.  The 
parameters  were  entered  into  the  Bidarsa  program  for  the  analysis.  The  resulting  bistatic 
maximum  range  product  is  74543  km  ,  and  the  estimated  RR  is  approximately  2  km. 
This  is  in  agreement  with  the  result  obtained  in  the  experiment  conducted. 


Satellite 


Target 


Direct 


Receiver  and  Processor 


Figure  5.2  NPS  bistatic  radar  experiment  configuration. 
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Effective  radiated  power  of  DBS- 1,  PtGt 

53  dBW 

Transmitting  frequency,/ 

12.7  GHz 

Effective  aperture  of  receiver  antenna,  Ae 

0.16  m^ 

Processor  gain  of  radar,  Gp 

61  dB 

Receiver  bandwidth.  Bn 

45  MHz 

Minimum  SNR  required  for  detection,  (S/N)min 

13  dB 

Distance  between  transmitter  and  receiver,  L 

37000  km 

Target  bistatic  RCS,  ob 

30  dBsm 

Two  way  Eosses,  E/r 

6dB 

Table  8.  System  parameters  for  the  NFS  bistatie  radar  experiment. 


The  experimental  result  confirmed  the  feasibility  of  the  concept,  and  suggested 
several  improvement  factors  [4]  that  could  employed  to  increase  the  detection  range. 
They  are: 

1 .  Using  of  low-noise  receiver  to  improve  the  SNR. 

2.  Increase  the  receiver  antenna  area  and  gain. 

3.  Increase  the  effective  radiated  power  on  the  transmitting  satellite. 

4.  Using  multiple  receivers  to  increase  the  processing  gain. 

It  is  important  to  take  note  that  in  the  bistatic-radar-satellite-hitchhiker  system,  each  6-dB 
increase  in  SNR  doubles  the  receiver-to-target  range,  Rr.  This  is  much  better  than  the 
monostatic  system,  whereby  it  takes  a  12-dB  increase  in  SNR  to  double  the  detection 
range. 


50 


D, 


SUMMARY 


In  this  chapter,  two  different  applieations  were  presented  and  the  Bidarsa  program 
was  used  to  carry  out  the  system  performance  prediction.  The  results  were  analyzed  and 
eompared  with  the  monostatic  radar  system.  Recommendations  were  also  diseussed  to 
improve  the  system  performance.  The  final  chapter  will  present  a  overall  summary  of 
this  thesis  and  suggest  possible  follow  up  work. 
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VI.  SUMMARY  AND  CONCLUSIONS 


A,  SUMMARY 

Bistatic  radar  offers  many  advantages  over  monostatic  radar  in  certain 
applications,  since  the  reeeiving  antenna  is  at  a  distant  site  with  reference  to  the 
transmitting  antenna.  This  thesis  achieved  the  objective  of  developing  a  software  model 
for  bistatic  radar  systems.  The  software  allows  a  user  to  examine  the  various  system 
design  parameters  and  evaluate  the  system  performance  tradeoffs,  under  different 
conditions  and  scenarios,  in  a  simulated  environment. 

MATLAB  GUI  is  chosen  as  the  software  to  develop  the  Bidarsa  simulation 
program  because  it  is  an  effective  interactive  tool  for  both  numerical  computation  and 
data  visualization.  It  also  allows  a  very  user-friendly  interface  to  be  designed  and 
operated. 

The  Bidarsa  program  is  easy  to  learn  and  navigate  with  Help  menus  designed  to 
guide  the  user  in  every  module.  It  is  suitable  to  be  used  by  any  engineering  student  or 
professionals  who  are  interested  in  the  bistatic  radar  field. 

With  the  eurrent  operation  trends  toward  a  lower  signature  target  and  the  third 
resurgence  for  bistatic  radar  using  Passive  Coherent  Location  [8],  the  bistatic  radar  offers 
a  much  greater  potential  now  than  at  any  time  before.  The  Bidarsa  program  will  help 
engineering  students  to  understand  the  fundamental  design  principles  of  bistatic  radar, 
and  provide  a  convenient  tool  for  the  professional  in  this  field  to  analyze  the  system 
performance.  The  software  is  available  free  of  eharge  and  ean  be  downloaded  from 
www.nps.navv.mil/ienn. 

B,  RECOMMENDATIONS 

While  this  thesis  focused  on  the  effort  to  develop  a  software  model  to  evaluate  the 
essential  performance  of  a  bistatic  radar  system,  there  are  other  areas  that  eould  be 
examined  to  provide  a  more  advanced  analysis  to  predict  the  bistatic  radar  performance. 
They  include; 
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1 .  Examine  the  effects  of  bistatic  RCS  for  surface  clutter,  which  consists  of 
both  ground  and  sea  echoes. 

2.  Clutter  tuning  concepts  that  cater  to  a  moving  transmitting  and  receiver. 

3.  Pulse  chasing  concepts  to  reduce  the  complexity  and  cost  of  multibeam 
bistatic  receivers. 

4.  Passive  Coherent  Location  (PCL)  considerations  [8],  such  as  the 
CELLDAR  concept  shown  in  Eigure  6.1.  CELLDAR  is  a  passive  system  based  on  the 
bistatic  radar  principle.  It  works  by  using  existing  transmissions  made  by  digital  cell 
phone  base  stations.  When  a  target  object  enters  the  detection  region,  these  cell  phone 
transmissions  are  reflected  by  it  and  detected  at  the  CELLAR  antennas. 


2.  Cell  phone 

signals  reflected 
by  aircraft  > 


1.  Aircraft  enters 
detection  region 


CELLDAR 

Concept  of  Operation 


Reflected  signals  & 
aircraft  noise  detected 
at  cellphone  radar 


Data  sent  to  Central  Control 
System  via  comms  network 


0  2002RoKje  Manor  Reaearch  Limited  I 


irii 


Cell  Phone 
Radar 


Acoustic 
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Sensor 
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Phone 

Base  Station 


Central 

Control  System 


Eigure  6.1  CELLDAR  concept  of  operation  proposed  by  Roke  Manor.  (Erom  Ref. 

[9].) 
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APPENDIX 


A,  MATLAB  CODE  FOR  BIDARSA.M 

This  appendix  details  the  Bidarsa.m  code  developed  to  allow  a  user  to  select  one 
of  the  three  group  functions:  (1)  specify  the  system  parameters,  (2)  carry  out  a  basic 
system  performance  analysis,  or  (3)  carry  out  an  extended  system  performance  analysis. 


function  varargout  =  Bidarsa(varargin) 

%  BIDARSA  M-file  for  Bidarsa.fig 

%  BIDARSA,  by  itself,  creates  a  new  BIDARSA  or  raises  the  existing 
%  singleton*. 

%  H  =  BIDARSA  returns  the  handle  to  a  new  BIDARSA  or  the  handle  to 
%  the  existing  singleton*. 

%  BIDARSA('CALLBACK',hObject,eventData,handles,...)  calls  the  local 
%  function  named  CALLBACK  in  BIDARSA.M  with  the  given  input  arguments. 
%  BIDARSA('Property','Value',.--)  creates  a  new  BIDARSA  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  Bidarsa  OpeningFunction  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 
%  stop.  All  inputs  are  passed  to  Bidarsa  OpeningFcn  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instance  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDEES 
%  Edit  the  above  text  to  modify  the  response  to  help  Bidarsa 
%  East  Modified  by  GUIDE  v2.5  20-0ct-2003  17:08:59 
%  Begin  initialization  code  -  DO  NOT  EDIT 
guiSingleton  =  I; 

gui_State  =  struct('gui_Name',  mfilename,  ... 

'gui_Singleton',  guiS  ingle  ton,  ... 

'gui_OpeningEcn',  @Bidarsa_OpeningEcn,  ... 

'gui_OutputEcn',  @Bidarsa_OutputEcn,  ... 

'gui  EayoutEcn',  []  ,  ... 

'gui_Callback',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callback  =  str2func(varargin{l}); 
end 

if  nargout 

[varargout) Lnargout}]  =  gui_mainfcn(gui_State,  varargin):}); 
else 

gui_mainfcn(gui_State,  varargin  ) : }  ); 
end 

%  End  initialization  code  -  DO  NOT  EDIT 
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%  —  Executes  just  before  Bidarsa  is  made  visible. 

function  Bidarsa_OpeningFcn(hObject,  eventdata,  handles,  varargin) 

handles,  output  =  hObject; 

guidata(hObject,  handles); 

clear  all; 

clc; 


%  —  Outputs  from  this  function  are  returned  to  the  command  line, 
function  varargout  =  Bidarsa_OutputFcn(hObject,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Executes  on  button  press  in  pushbuttonOO. 

function  pushbuttonOO_Callback(hObject,  eventdata,  handles) 

Defpara 

%  —  Executes  on  button  press  in  pushbutton  1 1 . 

function  pushbuttonl  l_Callback(hObject,  eventdata,  handles) 

RngCal 

%  —  Executes  on  button  press  in  pushbuttonl 2. 

function  pushbutton  12_Callback(hObject,  eventdata,  handles) 

TgtEoc 

%  —  Executes  on  button  press  in  pushbuttonl 3. 

function  pushbuttonl 3_Callback(hObject,  eventdata,  handles) 

CovAre 

%  —  Executes  on  button  press  in  pushbuttonl 4. 

function  pushbutton  14_Callback(hObject,  eventdata,  handles) 

FotPrt 

%  —  Executes  on  button  press  in  pushbutton21. 

function  pushbutton21_Callback(hObject,  eventdata,  handles) 

DopRel 

%  —  Executes  on  button  press  in  pushbutton3 1 . 

function  pushbutton31_Callback(hObject,  eventdata,  handles) 

EWEff 
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%  —  Executes  on  button  press  in  Close. 

function  Close_Callback(hObject,  eventdata,  handles) 

close  all  hidden 

%  —  Executes  on  button  press  in  Help. 

function  Help_Callback(hObject,  eventdata,  handles) 

yHlpMain 
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B, 


MATLAB  CODE  YOKDEFPARA.M 


This  appendix  details  the  Defpara.m  code  developed  to  provide  a  convenient  way 
for  user  to  specify  the  bistatic  radar’s  transmitter  and  receiver  design  parameters,  the 
target  information  and  the  environmental  information. 


function  varargout  =  Defpara(varargin) 

%  DEFPARA  M-file  for  Defpara.fig 

%  DEFPARA,  by  itself,  creates  a  new  DEFPARA  or  raises  the  existing 
%  singleton*. 

%  H  =  DEFPARA  returns  the  handle  to  a  new  DEFPARA  or  the  handle  to 
%  the  existing  singleton*. 

%  DEFPARA('CAEEBACK',hObject,eventData,handles,...)  calls  the  local 
%  function  named  CAEEBACK  in  DEFPARA.M  with  the  given  input  arguments. 
%  DEFPARA('Property','Value',.-.)  creates  a  new  DEFPARA  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  Defpara  OpeningFunction  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 
%  stop.  All  inputs  are  passed  to  Defpara  OpeningFcn  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instance  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDEES 
%  Edit  the  above  text  to  modify  the  response  to  help  Defpara 
%  East  Modified  by  GUIDE  v2.5  24-Oct-2003  01:25:10 
%  Begin  initialization  code  -  DO  NOT  EDIT 
gui_Singleton  =  I; 

gui_State  =  struct('gui_Name',  mfilename,  ... 

'guiSingleton',  guiS  ingle  ton,  ... 

'gui_OpeningFcn',  @Defpara_OpeningFcn,  ... 

'gui_OutputFcn',  @Defpara_OutputFcn,  ... 

'gui_EayoutFcn',  [] , ... 

'gui  Callback',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callback  =  str2func(varargin{l}); 
end 

if  nargout 

[varargout) Enargout}]  =  gui_mainfcn(gui_State,  varargin):}); 
else 

gui_mainfcn(gui_State,  varargin ) : }  ); 
end 

%  End  initialization  code  -  DO  NOT  EDIT 


%  —  Executes  just  before  Defpara  is  made  visible. 

function  Defpara_OpeningFcn(hObject,  eventdata,  handles,  varargin) 


58 


handles,  output  =  hObject; 
guidata(hObject,  handles); 
ele; 


%  —  Outputs  from  this  function  are  returned  to  the  command  line, 
function  varargout  =  Defpara_OutputFcn(hObject,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Pt_input_CreateFcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Pt_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Gt_input_CreateEcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Gt_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Ec[_input_CreateEcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Eq_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  SN_input_CreateEcn(hObject,  eventdata,  handles) 
if  ispc 


59 


set(hObject,'BackgroundColorVwhite'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  SN_input_Callback(hObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Gr_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Gr_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Ts_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Ts_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Bn_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Bn_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Br_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 
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set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Br_input_Callback(hObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Sm_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Sm_input_Callback(bObject,  eventdata,  bandies) 


function  Ba_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultFbcontrolBackgroundColor')); 

end 

function  Ba_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Ft_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultFbcontrolBackgroundColor')); 

end 

function  Ft_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Fr_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 
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function  Fr_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Lt_input_CreateFcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Lt_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Er_input_CreateEcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function  Er_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  on  button  press  in  Elpdate. 

function  Update_Callback(hObject,  eventdata,  handles) 

global  Pt  OO  Gt  OO  Eq_00  SN  OO  Pt  OO  Gr  OO  Ts  OO  Bn  OO  Br  OO... 

Sm  OO  Ba  OO  Ft  OO  Fr  OO  Et  OO  Er  OO  Eambda  OO 
clc; 

Pt_input  =  eval(get(handles.Pt_input,'String')); 

Gt_input  =  eval(get(handles.Gt_input,'String')); 

Fq_input  =  eval(get(handles.Fc[_input,'String')); 

SN_input  =  eval(get(handles.SN_input,'String')); 

Grinput  =  eval(get(handles.Gr_input,'String')); 

Ts_input  =  eval(get(handles.Ts_input,'String')); 

Bninput  =  eval(get(handles.Bn_input,'String')); 

Br_input  =  eval(get(handles.Br_input,'String')); 

Sminput  =  eval(get(handles.Sm_input,'String')); 

Ba_input  =  eval(get(handles.Ba_input,'String')); 

Ft_input  =  eval(get(handles.Ft_input,'String')); 

Fr_input  =  eval(get(handles.Fr_input,'String')); 

Etinput  =  eval(get(handles.Et_input,'String')); 

Er_input  =  eval(get(handles.Er_input,'String')); 
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Pt_00  =  10^(Pt_input/10); 

GtOO  =  10^(Gt_input/10); 

Fc[_00  =  Fq_input*le9; 

SN_00  =  (10^((SN_input-30)/10)); 

BaOO  =  deg2rad(Ba_input); 

GrOO  =  10^(Gr_input/10); 

TsOO  =  Tsinput; 

BnOO  =  Bn_input*le6; 

BrOO  =  Brinput; 

Sm_00  =  10^(Sm_input/10); 

Ft_00  =  10^(Ft_input/10); 

FrOO  =  10^(Fr_input/10); 

LtOO  =  10^(Lt_input/10); 

LrOO  =  10^(Lr_input/10); 

Lambda_00  =  3e8/(Fq_00);  %System  Wavelength 

mbO  =  msgbox(['  New  bistatie  radar  parameters  loaded.  '], 'Result'); 


%  —  Exeeutes  on  button  press  in  Print. 

function  Print_Callback(hObject,  eventdata,  handles) 

printdlg('-setup', handles. figurel) 

%  —  Executes  on  button  press  in  Close. 

function  Close_Callback(hObject,  eventdata,  handles) 

delete(handles .  figure  1 ) 

%  —  Executes  on  button  press  in  Help. 

function  Help_Callback(hObject,  eventdata,  handles) 

yHlpDefpara 
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C.  MATLAB  CODE  FOR  RNGCAL.M 

This  appendix  details  the  Rngcal.m  eode  developed  to  compute  the  bistatic 
maximum  range  product,  bistatic  radar  constant,  maximum  range-unambiguous  PRF, 
plotting  of  the  detection  contour,  and  carry  out  system  analysis  on  four  design 
parameters. 


function  varargout  =  RngCal(varargin) 

%  RNGCAL  M-file  for  RngCal.fig 

%  RNGCAL,  by  itself,  creates  a  new  RNGCAL  or  raises  the  existing 
%  singleton*. 

%  H  =  RNGCAL  returns  the  handle  to  a  new  RNGCAL  or  the  handle  to 
%  the  existing  singleton*. 

%  RNGCAL('CALLBACK',hObject,eventData,handles,...)  calls  the  local 
%  function  named  CALLBACK  in  RNGCAL.M  with  the  given  input  arguments. 
%  RNGCAL('Property','Value',.-.)  creates  a  new  RNGCAL  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  RngCal  OpeningFunction  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 
%  stop.  All  inputs  are  passed  to  RngCal  OpeningFcn  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instance  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDEES 
%  Edit  the  above  text  to  modify  the  response  to  help  RngCal 
%  East  Modified  by  GUIDE  v2.5  06-Nov-2003  08:56:10 
%  Begin  initialization  code  -  DO  NOT  EDIT 
gui_Singleton  =  I; 

gui_State  =  struct('gui_Name',  mfilename, ... 

'gui_Singleton',  guiS  ingle  ton,  ... 

'gui_OpeningEcn',  @RngCal_OpeningEcn,  ... 

'gui_OutputEcn',  @RngCal_OutputEcn,  ... 

'gui  EayoutEcn',  []  ,  ... 

'gui_Callback',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callback  =  str2func(varargin{l}); 
end 

if  nargout 

[varargout) Lnargout}]  =  gui_mainfcn(gui_State,  varargin):}); 
else 

gui_mainfcn(gui_State,  varargin  ) : }  ); 
end 

%  End  initialization  code  -  DO  NOT  EDIT 
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%  —  Executes  just  before  RngCal  is  made  visible. 

function  RngCal_OpeningFcn(hObject,  eventdata,  handles,  varargin) 

handles,  output  =  hObject; 

guidata(hObject,  handles); 

clc; 

axes(handles .  axes  1 ) ; 

[xl,map]  =  imread('z-Range.jpg',jpg'); 

image(xl); 

colonnap(map); 

set(handles.axesl,'VisibleVoff); 
global  PowP  FreP  RCSP  SNRP 
PowP  =  0; 

FreP  =  0; 

RCSP  =  0; 

SNRP  =  0; 


%  —  Outputs  from  this  function  are  returned  to  the  command  line, 
function  varargout  =  RngCal_OutputFcn(hObject,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Executes  on  mouse  press  over  figure  background,  over  a  disabled  or 

%  —  inactive  control,  or  over  an  axes  background. 

function  figure l_WindowButtonDownFcn(hObject,  eventdata,  handles) 


%  —  Executes  on  key  press  over  figure  1  with  no  controls  selected, 
function  figurel_KeyPressFcn(hObject,  eventdata,  handles) 


%  —  Executes  on  button  press  in  Calculate. 

function  Calculate_Callback(hObject,  eventdata,  handles) 

global  Pt  OO  Gt  OO  Fq_00  SN  OO  Pt  OO  Gr  OO  Ts  OO  Bn  OO  Br  OO... 

Sm  OO  Ba  OO  Ft  OO  Fr  OO  Et  OO  Er  OO  Eambda  OO  RtRrl  1  BRPl  1... 

PowP  FreP  RCSP  SNRP  ConP 
clc; 

k=  1.38065*10^(-23); 

Numll  =  (Pt_00*Gt_00*Gr_00*Eambda_00^2*Sm_00*Ft_00^2*Fr_00^2); 

Denll  =  ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Et_00*Er_00); 

RtRrl  1  =  ((Numl  1/Denl  1)^(0. 5))/1000;  %Bistatic  Max  Range  Product  in  km 

Denl2  =  ((4*pi)^3*k*Ts_00*Bn_00*Et_00*Er_00); 

BRPl  1  =  (Numl  1/Denl2);  %Bistatic  Radar  Parameter  in  km 

PRFull  =  3e8/((Br_00^2+(2*(RtRrll*1000)*(l+cos(Ba_00))))^0.5); 

PRFul2  =  PRFul  1/1000; 
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fprintf('\nBidarsa  -  Range  CalculationAn'); 

fprintf('Bistatic  Maximum  Range  Produet  RtRr  is  %-8.4e  km^2.\n',RtRrl  1); 
fprintf('Bistatic  Radar  eonstant  K  is  %-8.4e.\n\n',BRPl  1); 
fprintf('Maximum  Range-Unambiguous  PRF  is  %-8.4e  kHz.\n\n',PRFul2); 


close(figure(l)); 

%Plotting  Transmit  Power  Analysis 
if  (PowP==l) 

Ptll  =  (Pt_00/10):(Pt_00/10):(Pt_00M0); 

Numl  1 1  =  (Gt_00.*Pt_l  1  *Gr_00*Lambda_00^2*Sm_00*Ft_00^2*Fr_00^2); 

Benin  =  ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Lt_00*Lr_00); 

RtRrl  11  =  (((Numl  1 1  ./Deni  1 1).^(0.5))/1000);  %Bistatic  Max  Range  Product  in  km 

Pt_112  =  10.*logl0(Pt_ll); 

figure(l); 

plot(RtRrlll,Pt_112); 
title('Transmit  Power  Analysis'); 
xlabel('Max  Range  Product  (km^2)'); 
ylabel('Transmit  Power  (dBW)'); 
grid  on; 
end 

%Plotting  Tansmit  Frequency  Analysis 
if  (FreP==l) 

Fq_12  =  (Fq_00-5e9):(Fq_00/10):(Fq_00+5e9); 

Lambda_12  =  3e8./Fq_12; 

Numl21  =(Pt_00*Gt_00*Gr_00.*Lambda_12.^2*Sm_00*Ft_00^2*Fr_00^2); 
Denl21  =  ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Lt_00*Lr_00); 

RtRrl21  =  (((Numl21./Denl21).^(0.5))/1000);  %Bistatic  Max  Range  Product  in  km 

Fq_121  =Fq_12./le9; 

figure(l); 

plot(RtRrl21,Fq_121); 
title('Transmit  Frequency  Analysis'); 
xlabel('Max  Range  Product  (km^2)'); 
ylabel('Freq  (GHz)'); 
grid  on; 
end 


%Plotting  Target  RCS  Analysis 
if(RCSP==l) 

Sm_13  =  (Sm_00/10):(Sm_00/10):(Sm_00*10); 

Numl3 1  =  (Pt_00*Gt_00*Gr_00*Lambda_00^2.*Sm_13*Ft_00^2*Fr_00^2); 

Denl31  =  ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Lt_00*Lr_00); 

RtRrl31  =  (((Numl31./Denl31).^(0.5))/1000);  %Bistatic  Max  Range  Product  in  km 
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Sm_121  =  10.*logl0(Sm_13); 
figure(l); 

plot(RtRrl31,Sm_121); 
title('Target  RCS  Analysis'); 
xlabel('Max  Range  Product  (km^2)'); 
ylabel('Target  RCS  (dBsm)'); 
grid  on; 
end 

%Plotting  Receiver  SNR  Analysis 
if(SNRP==l) 

SN_14  =  (SN_00/10):(SN_00/10):(SN_00*10); 

Numl41  =  (Pt_00*Gt_00*Gr_00*Lambda_00^2*Sm_00*Ft_00^2*Fr_00^2); 

Denl41  =  ((4*pi)^3*k*Ts_00*Bn_00.*SN_14*Lt_00*Lr_00); 

RtRrl41  =  (((Numl41./Denl41)A(0.5))/1000);  %Bistatic  Max  Range  Product  in  km 
SN_141  =  (10.*logl0(SN_14))+30;  %Converting  dB  to  dBm 

figure(l); 

plot(RtRrl41,SN_141); 
title('SNR_{(min)}  Analysis'); 
xlabel('Max  Range  Product  (km^2)'); 
ylabel('SNR_{(min)}  (dBm)'); 
grid  on; 
end 


%  —  Executes  on  button  press  in  Pow_plot_input. 

function  Pow_plot_input_Callback(hObject,  eventdata,  handles) 


set(hObject,  'Value',  1); 
set(handles.Fre_plot_input,  'Value',  0); 
set(handles.RCS_plot_input,  'Value',  0); 
set(handles.SNR_plot_input,  'Value',  0); 
global  PowP  FreP  RCSP  SNRP 
PowP  =  1 ; 

FreP  =  0; 

RCSP  =  0; 

SNRP  =  0; 


%Turn  this  button  on 
%Turn  off  all  the  buttons 
%Turn  off  all  the  buttons 
%Turn  off  all  the  buttons 


%  —  Executes  on  button  press  in  Ere_plot_input. 
function  Ere_plot_input_Callback(hObject,  eventdata,  handles) 
set(hObject,  'Value',  1);  %Turn  this  button  on 

set(handles.Pow_plot_input,  'Value',  0);  %Turn  off  all  the  buttons 

set(handles.RCS_plot_input,  'Value',  0);  %Tum  off  all  the  buttons 

set(handles.SNR_plot_input,  'Value',  0);  %Tum  off  all  the  buttons 

global  PowP  EreP  RCSP  SNRP 
PowP  =  0;  %get(handles.Pow_plot_input,  'Value'); 
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FreP  =  1;  %get(hObject, 'Value'); 

RCSP  =  0; 

SNRP  =  0; 

%  —  Executes  on  button  press  in  RCS_plot_input. 
function  RCS_plot_input_Callback(bObject,  eventdata,  bandies) 
set(bObject,  'Value',  1);  %Turn  this  button  on 

set(bandles.Pow_plot_input,  'Value',  0);  %Turn  off  all  tbe  buttons 

set(bandles.Fre_plot_input,  'Value',  0);  %Tum  off  all  tbe  buttons 

set(handles.SNR_plot_input,  'Value',  0);  %Tum  off  all  the  buttons 

global  PowP  FreP  RCSP  SNRP 
PowP  =  0; 

FreP  =  0; 

RCSP  =  1; 

SNRP  =  0; 

%  —  Executes  on  button  press  in  SNR_plot_input. 
function  SNR_plot_input_Callback(hObject,  eventdata,  handles) 
set(hObject,  'Value',  1);  %Turn  this  button  on 

set(handles.Pow_plot_input,  'Value',  0);  %Tum  off  all  the  buttons 

set(handles.Fre_plot_input,  'Value',  0);  %Tum  off  all  the  buttons 

set(handles.RCS_plot_input,  'Value',  0);  %Tum  off  all  the  buttons 

global  PowP  FreP  RCSP  SNRP 
PowP  =  0; 

FreP  =  0; 

RCSP  =  0; 

SNRP=  1; 

%  —  Executes  on  button  press  in  Con_plot_input. 

function  Con_plot_input_Callback(hObject,  eventdata,  handles) 

close(figure(2)); 

CalculateCallback; 

Cass 


%  plot  ovals  of  cassini  for  bistatic  radar  SNR 

function  Cass 

global  BRPl  1  Br  OO 

rad=pi/180; 

i=0; 

F=Br_00; 

K=BRP11; 

delt=.l; 

a0=  10*logl0(16*K/(FM)); 
al  =  round(aO); 
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am6  =  al-6; 
am3  =  al-3; 
ap3  =  al+3; 
snr=[am6  am3  aO  ap3]; 

for  isnr=l;length(snr) 

SNR=  1 0^(snr(isnr)/ 1 0); 
i=i+l;  k=0; 

%  setup  for  b~=a 
a=L/2;  b4=K/SNR;  b=b4^(.25); 


if  b>=a 

for  thd=0:delt;360 
thr=thd*rad; 

%  b>a  is  one  curve 
if  b>a 
k=k+l; 

T(k)=thr; 

Rsq=a^2*(cos(2*thr)+sqrt(cos(2*thr)^2-l+(b/a)M)); 

rl(k)=sqrt(Rsq); 

r2(k)=0; 

end 

%  b=a  is  lemniscate 
if  b==a, 
k=k+l; 

T(k)=thr; 

Rsq=a^2*2*cos(2*thr);  rl(k)=sqrt(Rsq);  r2(k)=0; 
end 
end 

figure(2) 

if  i==l,  hold  on;  plot(-a,0,'or',a,0,'or');  end 
X=rl.*cos(T);  Y=rl.*sin(T); 
plot(X,Y); 
end 
end 

clear  rl  r2  T  X  Y 
%  b<a  is  two  curves 

%  this  covers  b<a  in  which  case  8  separate  curves  are  generated 
if  b<a 

for  curv=l;4 

thl=[0  90  180  270];  %lower  limit 

th2=[90  180  270  360];  %upper  limit 

k=0; 

%  decrease  angle  increment  for  these  small  curves 
for  thd=th  1  (curv) :  delt/ 1 0 ;  th2(curv) 
thr=thd*rad; 
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Rlsq=a^2*(cos(2*thr)+sqrt(cos(2*thr)^2-l+(b/a)^4)); 

R2sq=a^2*(cos(2*thr)-sqrt(cos(2*thr)^2-l+(b/a)M)); 

Rl=sqrt(Rlsq);  R2=sqrt(R2sq); 

if  imag(Rl)==0  &  imag(R2)==0  %Rl-'=  conj(R2) 

k=k+l; 

T(k)=tbr; 

rl(k)=(Rl);  r2(k)=(R2); 
end 
end 

figure(2);  X=rl.*eos(T);  Y=rl.*sin(T); 

plot(X,Y); 

elear  X  Y 

if  i==l,  bold  on;  end 
X=r2.*eos(T);  Y=r2.*sin(T); 
plot(X,Y); 
grid  on; 

elear  rl  r2  T  X  Y 
end 
end 

title('Deteetion  Contours  Analysis  -  Ovals  of  Cassini'); 

xlabel('Range  (km)'); 

ylabel('Range  (km)'); 

axesCposition',  [0  0  1  1], 'visible', 'off); 

text(0.31,0.48,'Tx'); 

text(0.70,0.48,'Rx'); 

bold  off 

fprintf('Ovals  of  Cassini  illustrate  tbe  eontour  plots  for  S/N  at:\n'); 
fprintf('%-4.2fdB,  %-4.2fdB,  %-4.2fdB  &  %-4.2fdB.\n\n',am6,  am3,  aO,  ap3); 


%  —  Exeeutes  on  button  press  in  Print. 

funetion  Print_Callbaek(bObjeet,  eventdata,  bandies) 

printdlgC-setup', bandies. figurel) 

%  —  Exeeutes  on  button  press  in  Close. 

funetion  Close_Callbaok(bObjeot,  eventdata,  bandies) 

delete(bandles .  figure  1 ) 

elose(figure(l)) 

elose(figure(2)) 

%  —  Exeeutes  on  button  press  in  Help. 

funetion  Help  Callbaek(bObieet,  eventdata,  bandies) 

yHlpRngCal 
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D,  MATLAB  CODE  FOR  TGTLOC.M 

This  appendix  details  the  TgtLoc.m  eode  developed  to  eompute  the  loeation  of  the 

target. 


funetion  varargout  =  TgtLoe(varargin) 

%  TGTLOC  M-file  for  TgtLoe.fig 

%  TGTLOC,  by  itself,  ereates  a  new  TGTLOC  or  raises  the  existing 
%  singleton*. 

%  H  =  TGTLOC  returns  the  handle  to  a  new  TGTLOC  or  the  handle  to 
%  the  existing  singleton*. 

%  TGTLOC('CALLBACK',hObjeot,eventData,handles,...)  ealls  the  loeal 
%  funetion  named  CALLBACK  in  TGTLOC.M  with  the  given  input  arguments. 
%  TGTLOC('Property','Value',...)  ereates  a  new  TGTLOC  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  TgtLoe  OpeningFunetion  gets  ealled.  An 
%  unreeognized  property  name  or  invalid  value  makes  property  applieation 
%  stop.  All  inputs  are  passed  to  TgtLoe  OpeningFen  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instanee  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDEES 
%  Edit  the  above  text  to  modify  the  response  to  help  TgtEoe 
%  East  Modified  by  GUIDE  v2.5  29-Oet-2003  15:06:14 
%  Begin  initialization  eode  -  DO  NOT  EDIT 
guiSingleton  =  1; 

gui_State  =  struot('gui_Name',  mfilename,  ... 

'guiSingleton',  guiS  ingle  ton,  ... 

'gui  OpeningEen',  @TgtEoo_OpeningEon,  ... 

'gui  OutputEen',  @TgtEoo_OutputEon,  ... 

'gui_EayoutEon',  []  ,  ... 

'gui_Callbaok',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callbaok  =  str2funo(varargin{l}); 
end 

if  nargout 

[varargout) Lnargout}]  =  gui_mainfon(gui_State,  varargin):}); 
else 

gui_mainfen(gui_State,  varargin  ) : }  ); 
end 

%  End  initialization  eode  -  DO  NOT  EDIT 


%  —  Exeeutes  just  before  TgtEoe  is  made  visible. 

funetion  TgtEoo_OpeningEon(hObjeot,  eventdata,  handles,  varargin) 

handles,  output  =  hObjeet; 
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guidata(hObject,  handles); 
global  TR  TT 
axes(handles .  axes  1 ) ; 

[xl,map]  =  imread('z-TgtLoe.jpg','jpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 

TR=  1; 

TT  =  0; 


%  —  Outputs  from  this  funetion  are  returned  to  the  eommand  line, 
funetion  varargout  =  TgtLoe_OutputFen(hObjeet,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  BR_input_CreateFon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion  BR_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  TR_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion  TR_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  TA_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 
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function  TA_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  on  button  press  in  Theta_R_input. 

function  Theta_R_input_Callback(hObject,  eventdata,  handles) 

set(hObject,  'Value',  1);  %Turn  this  button  on 

set(handles.Theta_T_input,  'Value',  0);  %Turn  off  all  the  buttons 

global  TR  TT 

TR  =  get(hObject,'Value'); 

TT  =  get(handles.Theta_T_input,  'Value'); 


%  —  Executes  on  button  press  in  Theta  T  input. 
function  Theta_T_input_Callback(hObject,  eventdata,  handles) 
set(hObject,  'Value',  1);  %Turn  this  button  on 

set(handles.Theta_R_input,  'Value',  0);  %Turn  off  all  the  buttons 

global  TR  TT 

TR  =  get(handles.Theta_R_input,  'Value'); 

TT  =  get(hObject,' Value'); 


%  —  Executes  on  button  press  in  Calculate. 

function  Calculate_Callback(hObject,  eventdata,  handles) 

global  Br  OO  TR  TT 

clc; 

TR_input  =  eval(get(handles.TR_input,'String')); 

TAinput  =  eval(get(handles.TA_input,'String')); 

BR_21  =  Br_00;  %Base  Range  in  km 

TR_21  =  TRinput; 

TA_21  =  deg2rad(TA_input); 

Elli  =BR_21/TR_21; 

Num21  =  BR_21*(1  -  E11C2); 

Den21  =  2*Elli*(l  +  Elli  *  sin(TA_21)); 

Rr2 1  =  Num2 1  /Den2 1 ;  %T gt  Range  to  Tx  in  km 

Num22  =  BR_21*(E11C2  +  1  +  (2*Elli*sin(TA_21))); 

Rt2 1  =  Num22/Den2 1 ;  %T gt  Range  to  Rx  in  km 

Den23  =  2*Elli*(l  -  Elli  *  sin(TA_21)); 

Rt22  =  Num21/Den23; 

Num24  =  BR_21*(E11C2  +  1  -  (2*Elli*sin(TA_21))); 

Rr22  =  Num24/Den23; 

fprintf('\nBidarsa  -  Target  EocationAn'); 
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ifTR>=TT 
Rt23  =  Rt21; 

Rr23  =  Rr21; 

Bet_21  =  asin(BR_21/Rt23*cos(TA_21)); 
fprintf('Theta_R  Selected\n'); 
else 

Rt23  =  Rt22; 

Rr23  =  Rr22; 

Bet_21  =  asin(BR_21/Rr23*cos(TA_21)); 
fprintf('Theta_T  Selected\n'); 
end 

fprintf('Distance  of  Target  to  Transmitter  is  %-8.4e  km.\n',Rt23); 
fprintf('Distanee  of  Target  to  Receiver  is  %-8.4e  km.\n',Rr23); 

Bet_22  =  rad2deg(Bet_21); 

fprintf('The  Beta  Angle  calculated  is  %-0.2f  deg.\n\n',  Bet_22); 


%  —  Executes  on  button  press  in  Print. 

function  Print_Callback(hObject,  eventdata,  handles) 

printdlg('-setup', handles. figurel) 

%  —  Executes  on  button  press  in  Close. 

function  Close_Callback(hObject,  eventdata,  handles) 

delete(handles .  figure  1 ) 

%  —  Executes  on  button  press  in  Help. 

function  Help_Callback(hObject,  eventdata,  handles) 

yHlpTgtEoc 
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E,  MATLAB  CODE  FOR  COVARE.M 

This  appendix  details  the  CovAre.m  code  developed  to  evaluate  the  coverage 
areas  of  both  detection-constrained  and  LOS-constrained  cases. 


function  varargout  =  CovAre(varargin) 

%  COVARE  M-file  for  CovAre.fig 

%  COVARE,  by  itself,  creates  a  new  COVARE  or  raises  the  existing 
%  singleton*. 

%  H  =  COVARE  returns  the  handle  to  a  new  COVARE  or  the  handle  to 
%  the  existing  singleton*. 

%  COVARE('CAEEBACK',hObject,eventData,handles,...)  calls  the  local 
%  function  named  CAEEBACK  in  COVARE.M  with  the  given  input  arguments. 
%  COVARE('Property','Value',...)  creates  a  new  COVARE  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  CovAre  OpeningFunction  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 
%  stop.  All  inputs  are  passed  to  CovAre  OpeningFcn  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instance  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDEES 
%  Edit  the  above  text  to  modify  the  response  to  help  CovAre 
%  East  Modified  by  GUIDE  v2.5  23-Oct-2003  1 1:49:58 
%  Begin  initialization  code  -  DO  NOT  EDIT 
guiSingleton  =  1; 

gui_State  =  struct('gui_Name',  mfilename,  ... 

'guiSingleton',  guiS  ingle  ton,  ... 

'gui_OpeningFcn',  @CovAre_OpeningFcn,  ... 

'gui_OutputFcn',  @CovAre_OutputFcn,  ... 

'gui_EayoutFcn',  []  ,  ... 

'gui_Callback',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callback  =  str2func(varargin{l}); 
end 

if  nargout 

[varargout) Enargout}]  =  gui_mainfcn(gui_State,  varargin):}); 
else 

gui_mainfcn(gui_State,  varargin  ) : }  ); 
end 

%  End  initialization  code  -  DO  NOT  EDIT 


%  —  Executes  just  before  CovAre  is  made  visible. 

function  CovAre_OpeningFcn(hObject,  eventdata,  handles,  varargin) 

handles,  output  =  hObject; 
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guidata(hObject,  handles); 
axes(handles .  axes  1 ) ; 

[xl,map]  =  imread('z-Centerr.jpgVjpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 

axes(handles .  axes2) ; 

[x2,map]  =  imread('z-Centert.jpgVjpg'); 

image(x2); 

eolormap(map); 

set(handles.axes2,'VisibleVoff); 

axes(handles.axes3); 

[x3,map]  =  imread('z-Centere.jpgVjpg'); 

image(x3); 

eolormap(map); 

set(handles.axes3,'VisibleVoff); 

axes(handles.axes4); 

[x4,map]  =  imread('z-LOS. jpg', ’jpg'); 

image(x4); 

eolormap(map); 

set(handles.axes4, 'Visible', ’off); 


%  —  Outputs  from  this  funetion  are  returned  to  the  eommand  line, 
funetion  varargout  =  CovAre_OutputFen(hObjeet,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Htg_input_CreateFon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor', 'white'); 
else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion  Htg_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Htx_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 
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else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

funetion  Htx_input_Callback(hObject,  eventdata,  bandies) 


%  —  Exeeutes  during  object  creation,  after  setting  all  properties, 
function  Hrx_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Hrx_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  on  button  press  in  Calculate. 

function  Calculate_Callback(bObject,  eventdata,  bandies) 

global  RtRrl  1  Br  OO 

clc; 

Htginput  =  eval(get(bandles.Htg_input,'String')); 
Htx_input  =  eval(get(bandles.Htx_input,'String')); 
Hrx_input  =  eval(get(bandles.Hrx_input,'String')); 


Htg_3 1  =  Htginput; 

Htx_3 1  =  Htx_input; 

Hrx_3 1  =  Hrxinput; 
BMRP31  =  RtRrl  1; 

BR31  =Br_00; 

TEST  =  2*sqrt(BMRP31); 
TEST2  =  abs(BR31-TEST); 
TEST3  =  TEST2/BR31; 


%Heigbt  in  km 


%Bistatic  Max  Range  Parameter  in  km 
%Base  Range  in  km 


fprmtf('\nBidarsa  -  Coverage  Area:\n'); 

if  (TEST3<=0.1)  %CheckifE~=2sqrt(k) 

Ab3 1  =  pi*BMRP3 1 ;  %Area  in  km^2 

fprintfCTarget  is  in  tbe  Eemniscate  Region.V); 
elseif(BR31<TEST) 

Ab31  =  pi*BMRP31*(l  -  (BR31M/BMRP3 1^2/64)  -  (3*BR31^8/BMRP31M/16384)); 
fprintfCTarget  is  in  tbe  Single  Cosite  Region.V); 
else  (BR31>TEST) 
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Ab31  =  (2*pi*BMRP31^2/BR31^2)*(l  +  (2*BMRP31^2/BR31M)  +  ... 

(12*BMRP31M/BR31^8)  +  (100*BMRP31^6/BR31^12)); 
fprintf('Target  is  in  the  Two  separate  Ellipses  Region.V); 
end 

fprintf('The  Deteetion-eonstrained  eoverage  area  is  %-8.4e  km^2.\n\n',Ab31); 


Rr3 1  =  130*(sqrt(Htg_3 1)  +  sqrt(Hrx_3 1));  %Range  in  km 
Rt31  =  130*(sqrt(Htg_31)  +  sqrt(Htx_31)); 

Theta_r  =  2*aeos((Rr31^2  -  Rt31^2  +  BR31^2)/(2*Rr31*BR31)); 
ThetaJ  =  2*aeos((Rt31^2  -  Rr31^2  +  BR31^2)/(2*Rt31*BR31)); 

if  (Rt31>=(BR31+Rr31)) 

Ae3 1  =  pi*Rr3 1^2;  %Area  in  km^2 

fprintf('Tx  eoverage  ineludes  all  Rx  eoverage  eirele.W); 
elseif  (Rr3 1  >=(BR3 1  +Rt3 1 )) 

Ae31  =  pi*Rt31^2; 

fprintf('Rx  eoverage  ineludes  all  Tx  eoverage  eirele.V); 
elseif  ((Rr3 1  +Rt3 1  )>=BR3 1 )  %Area  in  km^2 

Ae31  =  0.5*((Rr31^2*(Theta_r  -  sin(Theta_r)))  +  ... 

(Rt31^2*(Theta_t  -  sin(Theta_t)))); 
fprintf('Rx  and  Tx  is  in  LOS.V); 
else 

Ae31  =  0 

fprintf('No  LOS  between  Rx  and  Tx.V) 
end 

fprintf('The  LOS-eonstrained  coverage  area  is  %-8.4e  km^2.\n\n',Ac31); 


close(figure(l)); 

%Plotting  Transmitter-Target-Receiver  Height  Analysis 

BMRP32  =  RtRr  1 1  *  1  e3 ;  %Bistatic  Max  Range  Parameter  in  m 

BR32  =Br_00*le3;  %Base  Range  in  m 

for  L=  [(0.5*BR32)  BR32  (1.2*BR32)] 
n=0; 

for  ht=50;5;  1000  %target  heigth  in  m 

n=n+l; 

hrm=(  ((BMRP32+L^2/4)^0.5+L/2)/130-sqrt(ht*le3)  )^2; 
hr=hrm/1000; 

Ht(n)=ht; 

Hr(n)=hr; 

end 

figure(l) 

plot(Ht,Hr,'-'),  hold  on 
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title('Target  Height  Analysis'); 
xlabel('Target  Altitude  (m)') 
ylabel('Tx  or  Rx  Antenna  Altitude  (m)') 
grid  on 
end 

axesCposition',  [0  0  1  1], 'visible', ’off); 

text(0.30,0.15,'0.5*L'); 

text(0.48,0.20,'L'); 

text(0.68,0.26,'1.2*L'); 

hold  off 


%  —  Executes  on  button  press  in  Print. 

function  Print_Callback(hObject,  eventdata,  handles) 

printdlg('-setup', handles. figurel) 

%  —  Executes  on  button  press  in  Close. 

function  Close_Callback(hObject,  eventdata,  handles) 

delete(handles .  figure  1 ) 

close(figure(l)) 

%  —  Executes  on  button  press  in  Help. 

function  Help_Callback(hObject,  eventdata,  handles) 

yHlpCovAre 
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F.  MATLAB  CODE  FOR  FOTPRT.M 

This  appendix  details  the  ForPrt.m  eode  developed  to  evaluate  the  bistatie 
footprint  and  elutter  area. 


funetion  varargout  =  FotPrt(varargin) 

%  FOTPRT  M-file  for  FotPrt.fig 

%  FOTPRT,  by  itself,  ereates  a  new  FOTPRT  or  raises  the  existing 
%  singleton*. 

%  H  =  FOTPRT  returns  the  handle  to  a  new  FOTPRT  or  the  handle  to 
%  the  existing  singleton*. 

%  FOTPRT('CALLBACK',hObjeot,eventData,handles,...)  ealls  the  loeal 
%  funetion  named  CALLBACK  in  FOTPRT.M  with  the  given  input  arguments. 
%  FOTPRT('Property','Value',...)  ereates  a  new  FOTPRT  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  FotPrt  OpeningFunetion  gets  ealled.  An 
%  unreeognized  property  name  or  invalid  value  makes  property  applieation 
%  stop.  All  inputs  are  passed  to  FotPrt  OpeningFen  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instanee  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDEES 
%  Edit  the  above  text  to  modify  the  response  to  help  EotPrt 
%  East  Modified  by  GUIDE  v2.5  20-0et-2003  17:10:38 
%  Begin  initialization  eode  -  DO  NOT  EDIT 
guiSingleton  =  1; 

gui_State  =  struot('gui_Name',  mfilename,  ... 

'guiSingleton',  guiS  ingle  ton,  ... 

'gui_OpeningEon',  @EotPrt_OpeningEon,  ... 

'gui  OutputEen',  @EotPrt_OutputEon,  ... 

'gui_EayoutEon',  []  ,  ... 

'gui_Callbaok',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callbaok  =  str2funo(varargin{l}); 
end 

if  nargout 

[varargout) Lnargout}]  =  gui_mainfon(gui_State,  varargin):}); 
else 

gui_mainfen(gui_State,  varargin  ) : }  ); 
end 

%  End  initialization  eode  -  DO  NOT  EDIT 


%  —  Exeeutes  just  before  EotPrt  is  made  visible. 

funetion  PotPrt_OpeningPon(hObjeot,  eventdata,  handles,  varargin) 

handles,  output  =  hObjeet; 
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guidata(hObject,  handles); 
axes(handles .  axes  1 ) ; 

[x,map]  =  imread('z-BWLimit.jpgVjpg'); 

image(x); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 

axes(handles .  axes2) ; 

[y,map]  =  imread('z-PWLimit.jpgVjpg'); 

image(y); 

eolormap(map); 

set(handles.axes2,'VisibleVoff); 


%  —  Outputs  from  this  funetion  are  returned  to  the  eommand  line, 
funetion  varargout  =  FotPrt_OutputFen(hObjeet,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Htx_input_CreateFon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion  Htx_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Hrx_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion  Hrx_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Rr_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 


81 


set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Rr_input_Callback(hObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Rt_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Rt_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Rx_BW_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Rx_BW_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Tx_BW_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Tx_BW_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  PusWid_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 
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function  PusWid_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  on  button  press  in  Calculate. 

funetion  Calculate_Callbaek(hObject,  eventdata,  handles) 

global  Ba  OO  Br  OO 

ele; 

Rrxinput  =  eval(get(handles.Rr_input, 'String')); 

Rtxinput  =  eval(get(handles.Rt_input, 'String')); 

RBWinput  =  eval(get(handles.Rx_BW_input, 'String')); 

TBWinput  =  eval(get(handles.Tx_BW_input, 'String')); 

PuW_input  =  eval(get(handles.PusWid_input, 'String')); 

Bet_41  =  Ba_00;  %Beta  Angle  in  rad 

Rrx_41  =  Rrxinput; 

Rtx_41  =  Rtxinput; 

RBW_41  =  deg2rad(RBW_input); 

TBW_41  =  deg2rad(TBW_input); 

PuW_41  =  PuW_input*  10^-6; 

Num4 1  =  Rrx_4 1  *RBW_4 1  *Rtx_4 1  *TBW_4 1 ; 

Den41  =  sin(Bet_41); 

Num42  =  3e8*PuW_41*Rrx_41*RBW_41; 

Den42  =  2*(eos(Bet_41/2))^2; 

DelR  =  (3e8*PuW_41)/(2*cos(Bet_41/2)); 

RrDr  =  Rrx_41*RBW_41; 

CeS41  =  3e8*PuW_41/(2*eos(Ba_00/2)); 

Br41  =Br_00*le3; 

Aa41  =  (Rrx_41+Rtx_41)*le3/2; 

Ap41  =  Aa41  +  (3e8*PuW_41/2); 

Bb41  =  (Aa41^2  -  (Br41^2)/4)^(0.5); 

Bp41  =  (Ap41^2  -  (Br41^2)/4)^(0.5); 

Emax41  =  ((Aa41*(Ap41-Aa41))/(Bb41*(Bp41-Bb41 )))-!; 

fprintf('\nBidarsa  -  Bistatic  EootprintAn'); 
if  DelR  <  RrDr 
Ae41  =  Num42/Den42; 
fprintf('This  is  a  Pulsewidth  Eimited  case.\n'); 
else 
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Ac41  =  Num41/Den41; 
fprintf('This  is  a  Beamwidth  Limited  caseAn'); 
end 

fprintf('The  Bistatie  Footprint  is  %-8.4e  sqr-kmAn\n',Ao41); 

fprintf('Range  Cell  Size  is  %-8.4e  m.\n',CeS41); 
fprintf('Maximum  Error  is  %-8.4e  mAn\n',Emax41); 


%  —  Executes  on  button  press  in  Print. 

function  Print_Callback(hObject,  eventdata,  handles) 

printdlg('-setup', handles. figurel) 

%  —  Executes  on  button  press  in  Close. 

function  Close_Callback(hObject,  eventdata,  handles) 

delete(handles .  figure  1 ) 

%  —  Executes  on  button  press  in  Help. 

function  Help_Callback(hObject,  eventdata,  handles) 

yHlpEotPrt 
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G.  MATLAB  CODE  FOR  DOPREL.M 

This  appendix  details  the  DopRel.m  code  developed  to  examine  the  doppler  shift 
of  a  moving  target. 


function  varargout  =  DopRel(varargin) 

%  DOPREL  M-fde  for  DopRel.fig 

%  DOPREL,  by  itself,  creates  a  new  DOPREL  or  raises  the  existing 
%  singleton*. 

%  H  =  DOPREL  returns  the  handle  to  a  new  DOPREL  or  the  handle  to 
%  the  existing  singleton*. 

%  DOPREL('CALLBACK',hObject,eventData,handles,...)  calls  the  local 
%  function  named  CALLBACK  in  DOPREL.M  with  the  given  input  arguments. 
%  DOPREL('Property','Value',.--)  creates  a  new  DOPREL  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  DopRel  OpeningLunction  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 
%  stop.  All  inputs  are  passed  to  DopRel  OpeningLcn  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instance  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDLES 
%  Edit  the  above  text  to  modify  the  response  to  help  DopRel 
%  Last  Modified  by  GUIDE  v2.5  02-Nov-2003  11:11:15 
%  Begin  initialization  code  -  DO  NOT  EDIT 
gui_Singleton  =  I; 

gui_State  =  struct('gui_Name',  mfilename,  ... 

'guiSingleton',  guiS  ingle  ton,  ... 

'gui  OpeningEcn',  @DopRel_OpeningEcn,  ... 

'gui  OutputEcn',  @DopRel_OutputEcn,  ... 

'gui  LayoutEcn',  []  ,  ... 

'gui  Callback',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callback  =  str2func(varargin{l}); 
end 

if  nargout 

[varargout) Enargout}]  =  gui_mainfcn(gui_State,  varargin):}); 
else 

gui_mainfcn(gui_State,  varargin ) : }  ); 
end 

%  End  initialization  code  -  DO  NOT  EDIT 


%  —  Executes  just  before  DopRel  is  made  visible. 

function  DopRel_OpeningEcn(hObject,  eventdata,  handles,  varargin) 

handles,  output  =  hObject; 
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guidata(hObject,  handles); 
axes(handles .  axes  1 ) ; 

[xl,map]  =  imread('z-DopRel.jpgVjpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 


%  —  Outputs  from  this  funetion  are  returned  to  the  eommand  line, 
funetion  varargout  =  DopRel_OutputFen(hObjeet,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Vel_input_CreateFon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion  Vel_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Del_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion  Del_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  TA_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion  TA_input_Callbaok(hObjeot,  eventdata,  handles) 
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%  —  Executes  on  button  press  in  Calculate. 

function  Calculate_Callback(hObject,  eventdata,  handles) 

global  Ba  OO  Lambda  OO 

clc; 

Velinput  =  eval(get(handles.Vel_input, 'String')); 

Delinput  =  eval(get(handles.Del_input, 'String')); 

Vel_51  =  Velinput; 

Del_5 1  =  deg2rad(Del_input); 

DoS_51  =  (2*Vel_51/Lambda_00)*(cos(Del_51)*cos(Ba_00/2)); 
DoS_52  =  DoS_5 1/1000;  %in  kHz 

fprintf('\nBidarsa  -  Doppler  RelationAn'); 
fprintf('Bistatic  Doppler  Shift  is  %-8.4e  kHz.\n\n',DoS_52); 


close(figure(l)); 

%Plotting  Doppler  Shift  Analysis 
forBa_51  =  0:(pi/4):pi; 

Del_52  =  -pi:(pi/360):pi; 

DoS_53  =  (2*Vel_51/Lambda_00).*cos(Del_52).*cos(Ba_51/2); 
figure(l); 

DoS_54  =  DoS_53./le3; 

Del_53  =  l.*rad2deg(Del_52); 
plot(Del_53,DoS_54); 

title('Bistatic  Target  Doppler  Shift  Analysis'); 
xlabel('Target  velocity  aspect  angle,  \delta,  (deg)'); 
ylabel('Bistatic  Doppler  Shift,  F_{tgt}(kHz)'); 
xhm([-180  180]); 
grid  on; 
hold  on 
end 

axesCposition',  [0  0  1  l],'visible','off); 
text(0.87,0.54,'\beta  =  180\circ'); 
text(0.87,0.40,'\beta  =  135\circ'); 
text(0.87,0.25,'\beta  =  90\circ'); 
text(0.87,0.17,'\beta  =  45\circ'); 
text(0.87,0.13,'\beta  =  0\circ'); 
hold  off 


%  —  Executes  on  button  press  in  Print. 

function  Print_Callback(hObject,  eventdata,  handles) 

printdlg('-setup',handles.figurel) 
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%  —  Executes  on  button  press  in  Close. 

function  Close_Callback(hObject,  eventdata,  handles) 

delete(handles .  figure  1 ) 

close(figure(l)) 

%  —  Executes  on  button  press  in  Help. 

function  Help_Callback(hObject,  eventdata,  handles) 

yHlpDopRel 
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H.  MATLAB  CODE  FOR  EWEFF.M 

This  appendix  details  the  EWEjf.m  code  developed  to  evaluate  the  bistatic  radar 
performance  under  different  jammer  environments. 


function  varargout  =  EWEff(varargin) 

%  EWEEE  M-fde  for  EWEff  fig 

%  EWEEE,  by  itself,  creates  a  new  EWEEE  or  raises  the  existing 
%  singleton*. 

%  H  =  EWEEE  returns  the  handle  to  a  new  EWEEE  or  the  handle  to 
%  the  existing  singleton*. 

%  EWEEE('CAEEBACK',hObject,eventData,handles,...)  calls  the  local 
%  function  named  CAEEBACK  in  EWEEE. M  with  the  given  input  arguments. 
%  EWEEE('Property','Value',.--)  creates  a  new  EWEEE  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 
%  applied  to  the  GUI  before  EWEff  OpeningEunction  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 
%  stop.  All  inputs  are  passed  to  EWEff  OpeningEcn  via  varargin. 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instance  to  run  (singleton)". 

%  See  also:  GUIDE,  GUIDATA,  GUIHANDEES 
%  Edit  the  above  text  to  modify  the  response  to  help  EWEff 
%  East  Modified  by  GUIDE  v2.5  23-Oct-2003  14:10:38 
%  Begin  initialization  code  -  DO  NOT  EDIT 
gui_Singleton  =  I; 

gui_State  =  struct('gui_Name',  mfilename,  ... 

'guiSingleton',  guiS  ingle  ton,  ... 

'gui_OpeningEcn',  @EWEff_OpeningEcn,  ... 

'gui_OutputEcn',  @EWEff_OutputEcn,  ... 

'gui  EayoutEcn',  []  ,  ... 

'gui  Callback',  []); 
if  nargin  &  isstr(varargin{l}) 

gui_State.gui_Callback  =  str2func(varargin{l}); 
end 

if  nargout 

[varargout) Enargout}]  =  gui_mainfcn(gui_State,  varargin):}); 
else 

gui_mainfcn(gui_State,  varargin ) : }  ); 
end 

%  End  initialization  code  -  DO  NOT  EDIT 


%  —  Executes  just  before  EWEff  is  made  visible. 

function  EWEff_OpeningEcn(hObject,  eventdata,  handles,  varargin) 

handles,  output  =  hObject; 
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guidata(hObject,  handles); 
axes(handles .  axes  1 ) ; 

[xl,map]  =  imread('z-EWEff.jpgVjpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 


%  —  Outputs  from  this  funetion  are  returned  to  the  eommand  line, 
funetion  varargout  =  EWEff_OutputEen(hObjeet,  eventdata,  handles) 
varargoutjl}  =  handles,  output; 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Pj_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEheontrolBaekgroundColor')); 

end 

funetion  Pj_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Gj_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion  Gj_input_Callbaok(hObjeot,  eventdata,  handles) 


%  —  Exeeutes  during  objeet  ereation,  after  setting  all  properties, 
funetion  Ej_input_CreateEon(hObjeot,  eventdata,  handles) 
if  ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion  Ej_input_Callbaok(hObjeot,  eventdata,  handles) 
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%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Bn_input_CreateFcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColorVwhite'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Bn_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Ej_input_CreateEcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function  Ej_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Rj_input_CreateEcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function  Rj_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Bj_input_CreateEcn(hObject,  eventdata,  handles) 
if  ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function  Bj_input_Callback(hObject,  eventdata,  handles) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Sm_input_CreateEcn(hObject,  eventdata,  handles) 
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if  ispc 

set(hObject,'BackgroundColorVwhite'); 

else 

set(hObjeet,'BackgroundColor',get(0,'defaultUicontrolBaekgroundColor')); 

end 

function  Sm_input_Callback(hObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Lr_input_CreateFcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function  Lr_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Et_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Et_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Er_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function  Er_input_Callback(bObject,  eventdata,  bandies) 


%  —  Executes  during  object  creation,  after  setting  all  properties, 
function  Et_input_CreateEcn(bObject,  eventdata,  bandies) 
if  ispc 

set(bObject,'BackgroundColor','wbite'); 
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else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

funetion  Ft_input_Callbaek(hObjeet,  eventdata,  bandies) 


%  —  Executes  on  button  press  in  Calculate. 

function  Calculate_Callback(bObject,  eventdata,  bandies) 

global  Pt  OO  Gt  OO  Fq_00  SN  OO  Pt  OO  Gr  OO  Ts  OO  Bn  OO  Sm  OO... 

Ft  OO  Fr  OO  Ft  OO  Fr  OO  Fambda  OO  RtRrl  1 
clc; 

Pj_input  =  eval(get(bandles.Pj_input, 'String')); 

Gj_input  =  eval(get(bandles.Gj_input, 'String')); 

Fjinput  =  eval(get(bandles.Fj_input, 'String')); 

Bj_input  =  eval(get(bandles.Bj_input, 'String')); 

Rj_input  =  eval(get(bandles.Rj_input, 'String')); 

Fj_input  =  eval(get(bandles.Fj_input, 'String')); 

Pj_61  =  10^(Pj_input/10); 

Gj_61  =  10^(Gj_input/10); 

Fj_61  =  10^(Fj_input/10); 

Bj_61  =  Bj_input*10^6; 

Rj_61  =  Rj_input*10^3; 

Fj_61  =  10^(Fj_input/10); 

Num61  =  Pt_00*Gt_00*(Ft_00)^2*(Fr_00)^2*Fj_61  *Bj_61  *Sm_00*(Rj_61)^2; 

Dem61  =  Pj_61*Gj_61*(Fj_61)^2*Ft_00*Fr_00*Bn_00*4*pi*SN_00; 

RtRrJl  =  (Num61/Dem61)^(0.5); 
fprintf('\nBidarsa  -  EW  Effects  :\n'); 

fprintf('Bistatic  Max  Range  Product  with  Jamming  (RtRr)j  is  %-8.4e  km^2.\n\n', RtRrJl); 


close(figure(l)); 

%Plottmg  EW  Effects  Analysis 

for  Sm_611  =  [(Sm  OO/lO)  Sm  OO  (Sm_00*10)] 

Pj_611  =  (Pj_6 1/1 00):(Pj_6 1/1 00):(Pj_6 1*100); 

Num61 1  =  Pt_00*Gt_00*(Ft_00)^2*(Fr_00)^2*Ej_61  *Bj_61  .*Sm_61 1  *(Rj_61)^2; 
Dem61 1  =  Gj_61.*Pj_61  l*(Fj_61)^2*Et_00*Er_00*Bn_00*4*pi*SN_00; 

RtRrJ61  =  ((Num61 1./Dem61 1).^(0.5)); 

Pj_612  =  10.*logl0(Pj_611); 
figure(l); 

plot(RtRrJ61,Pj_612); 
title('Burn-tbrougb  Range  Analysis'); 


93 


xlabel('Max  Range  Product  with  Jamming  (km^2)'); 
ylabel('Jammer  Transmit  Power  (dBW)'); 
grid  on; 
hold  on 
end 

axesCposition',  [0  0  1  l],'visibleVoff); 
text(0 .15,0. 15, '0.1  \sigma_  {B  } '); 
text(0.38,0.20,'\sigma_{B}'); 
text(0.60,0.25,'10\sigma_{B}'); 
hold  off 


%  —  Executes  on  button  press  in  Print. 

function  Print_Callback(hObject,  eventdata,  handles) 

printdlg('-setup',handles.figurel) 

%  —  Executes  on  button  press  in  Close. 

function  Close_Callback(hObject,  eventdata,  handles) 

delete(handles .  figure  1 ) 

close(figure(l)) 

%  —  Executes  on  button  press  in  Help. 

function  Help_Callback(hObject,  eventdata,  handles) 

yHlpEWEff 
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